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1 Introduction 

Discovered in 1982 by Tsui, Stormer and Gossard the fractional Quantum Hall effect 
(FQHE) has opened a field of vivid activity until today. The reason for its continuing actuality 
can be attributed mainly to four fascinating features of this two-dimensional electron gas in 
a magnetic field. 

As in the integer quantum Hall effect (IQHE) |2l the measured Hall resistance in the 
FQHE is reproduced to very high accuracy due to a topological origin of the quantization. This 
was already pointed out by Laughlin [2j short after the discovery of the effect. Quantization 
that can be traced back to topological features is not restricted to the IQHE and FQHE but 
was also found in superfluidity, superconductivity and in Josephson junctions. A collection 
on this subject by Thouless can be found in [5]. Secondly, the FQHE is an effect caused by 
correlations. It can be described by a Hamiltonian solely containing two-particle interactions 
(see section ll.l.4|) . Thus, it belongs to the most strongly correlated systems studied so far. 
Closely linked to the correlations is the appearance of fractionally charged quasiparticles first 
introduced by Laughlin The prediction of these exotic particles inspired the efforts of 
experimentalists and theorists and the particles' existence became more and more manifest. 
Nevertheless investigating their properties - such as in tunneling experiments - is still a do- 
main of recent experiments jl9l l2Uj and accompanying theoretical work jl7j . These tunneling 
experiments raised the question what the dynamics of quasiparticles at a tunneling constric- 
tion is. It is also the motivation for this diploma thesis. The last point is a more theoretical 
one: In some cases the edges of a sample in the FQHE ~ like in the IQHE - can be considered 
to be an embodiment of an one-dimensional interacting system for which analytically solvable 
models [Luttinger liquid) are already known independently. 

Motivated by experiments El HO] and theory ^1 El on single quasiparticle tunneling 
through a quantum point contact, a constricted fractional quantum Hall system will also 
be investigated in this work (section |H1) . To address the question of a single quasiparticle 
tunneling through this constriction, here we resort to quasiholes for which well approved trial 
wavefunctions are known. Basing the work on the electronic many-particle Hamiltonian and 
these trial wavefunctions allows for a view on quasihole tunneling that is independent on the 
chiral Luttinger liquid model used so far for explaining the experiments. Finite systems of 
few (4 - 6) electrons in the FQH regime will be investigated, which permits the electronic 
many-particle Hamiltionian to be dealt with by exact numerical diagonalization. 

The outline of the main part of this work, which is divided into six sections, is as follows. 
Section 12 covers the major developments on the fractional quantum Hall regime concerning 
this work. In doing so, the question of quasiparticle tunneling will be touched and embedded 
into past and ongoing research. In section |21 the single particle basis for the chosen boundary 
conditions will be derived. The connection of the boundary conditions to an electric in-plane 
field will also be pointed out. Section covers the homogeneous (absent external potential) 
many-electron system. A short-ranged interaction will be introduced and shown to be useful 
for the following work. The features of a system at a filling factor v = ^ will be reviewed 
for Coulomb interaction and this short-range interaction. Two different methods of inserting 
a quasihole into the system will be derived from Laughlin's trial wavefunctions |21 and the 
stability of these excitations will be checked for both electron-electron interactions. Aside 
from that, a system providing bound quasihole states is found. These bound states will 
be used in section to create the most simple system in which tunneling of quasiholes can 
be observed. In section j^l an external potential is introduced to model a quantum point 
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contact. Corrections to the current operators arising from contributions of the next Landau 
level turn out to be crucial to obtain a consistent picture. The time evolution of an injected 
quasihole will be evaluated for both, weak and strong external potentials (compared to the 
excitation gap). Creating a tunneling constriction by a strong potential is counteracted by 
the incompressibility of the system. Ways to overcome this competition are examined and 
lead to a realization of an effective tunneling barrier in which the time evolution of a quasihole 
can be studied. The last section contains a summary of the main results and conclusions. 
Perspectives on possible further investigations and on questions not exhausted in this work 
are given as well. 
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2 Quasiparticles in past and recent research: An overview on tiie 
FQH regime 

In the following, the major developments on the fractional quantum Hall regime from the 
perspective of this work will be mentioned. This is not only intended to give a coarse overview 
but also to render more precise the question about quasiparticle tunneling and to relate it to 
past and recent research. 

As stated above, it was Laughlin [HI who found the low energy excitations in a FQHE 
system to carry a fraction of an electron's charge. The finite creation energy of these quasi- 
particles explained the crucial physics like the incompressibility (section I4.1.5j) which in turn 
helped explaining the vanishing longitudinal resistance px^x = 0. His trial wavefunction also 
identified the origin of the gapped ground state (section 14. 1|) and revealed the uniqueness 
of the filling factors v = m odd. Moreover, the one-to-one correspondence between flux 
quanta and zeros in the many-body wavefunctions (vortices) became apparent. Numerical 
work by Yoshioka, Halperin and Lee on finite systems with rectangular geometry was in 
agreement with Laughlin's proposal and also finite systems with spherical symmetry investi- 
gated by Haldane and Rezayi jZj confirmed Laughlin's wavefunction for the ground state and 
the quasiparticle excitations. 

The filling factors v = 2pm+i PjItt- ^ where the FQHE was observed as well, could not be 
explained by Laughlin's theory. Haldane explained it in a hierachichal scheme [HI as the FQHE 
of quasiparticles. Again the quasiparticles played an essential role. Another finding in this 
paper was the importance of the short-ranged part of the Coulomb interaction that was proven 
to be responsible for the appearance of the correlations in Laughlin's wavefunction. This was 
precised by Trugman and Kievelson [21] who even showed Laughlin's wavefunction to be an 
exact eigenstate for a special kind of short-ranged electron-electron interaction (section 14. Ij) . 

These earlier developments which base directly on the many-particle Hamiltonian of the 
system (section 13. Ij) together with the trial wavefunctions of fractionally charged excitations 
fsection l4.2() constitute the footing for the work at hand. This is to point out the independence 
of our results on other theoretical work resting upon more elaborate theoretical models. 
Although not used explicitly in this work, these developments shall be sketched here because 
they are mainly used in explaining the most recent experiments of quasiparticle tunneling. 

Jain introduced a new quasiparticle - a compound of an electron and an even number of 
vortices ~ the composite fermion jH] . This invention made a unification of the IQHE and FQHE 
possible: "The FQHE is the IQHE of composite fermions". This unification also includes an 
alternative to the hierachical picture in explaining the fractions v = 2pm+i • According to 
Jain, the correlations are the key issue of this concept: The interaction between the electrons 
is incorporated into the system in the definition of the composite fermions. The system of 
strongly correlated electrons thus transforms into a system of weakly interacting composite 
fermions. Their fermionic nature attended by weak interactions among them allow these 
compound particles to show the IQHE again. Due to the weak interaction, this invention 
opened the field towards the well approved mean field descriptions and those that go beyond. 
An overview on the wealth of developments basing on composite fermions can be found in [Hj. 
The success of this theory proved the importance of two kinds of correlations for the FQHE: 
The Laughlin-like correlations that bind vortices to electrons resident in the definition of the 
composite fermions and the fermionic Pauli correlations responsible for the rigidity of effect. 

Recent experiments on quasiparticle tunneling |18[ ll9 | l2fl] are closely related to a description 
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of the fractional quantum Hah regime by edge states. This theory, used successfuUy in 
explaining the IQHE, was ported to the FQH regime. Its applicability in the latter case 
can intuitively be understood due to the unification of IQHE and FQHE by the analogy of 
electrons and composite fermions. 

A more rigor treatment is based on the wavefunction picture. In an abruptly confined two- 
dimensional electron gas at integer filling the only excitations allowed by Pauli's principle 
appear at the edges of the sample where the Landau level crosses the Fermi energy. In the 
wavefunction picture for the IQHE these excitations can be shown to be bosonic. The same 
is true for the FQHE JU]. Creating an effective low energy theory of excitations near the 
two Fermi points of the one-dimensional edge makes the Luttinger liquid model applicable in 
the case of the IQHE. Similar low energy bosonic excitations in the FQHE are treated in the 
chiral Luttinger liquid model developed by Wen . A subsequent work of Chamon and Wen 
|12j pointed out the importance of the confining potential to produce abrupt edges. "Soft" 
potentials were shown to produce more complicated edge structures. 

Quasiparticle tunneling into the edges is considered to be a proper means of verifying the 
chiral Luttinger liquid nature and of investigating the structure of the edges. Such tunneling 
experiments, first used for examining the quasiparticles themselves rather than the edges, were 
performed by Simmons et. al. J3] already before Wen's theory. They measured fluctuations 
in the longitudinal conductivity pxx which were believed to origin from tunneling processes of 
quasiparticles between the current carrying edges as explained in a paper by Pokrovsky and 
Pryadko However, the fractional filling of the system could not be ruled out completely 
as a source of these fluctuations. To exclude them, non-equilibrium processes had to be 
investigated. 

Kane and Fisher ^S] developed a non-equilibrium theory that allowed the calculation of 
fluctuations in the current through a FQH system. It is based on the chiral Luttinger liquid 
model by Wen and describes the tunneling process as a point-like tunneling impurity 
between two edges of the sample. To sketch the idea, a schematics of the two edge channels 
that conduct the current through the sample is given in Fig. ^ These edge channels are bent 
due to an external potential V{x, y) that creates a quantum point contact in the system. The 
quasiparticles flowing along the edges are believed to tunnel from one edge to the other and 
therefore cause fluctuations in time in the current /. These can be directly observed in the 
experiment. 

An experiment of a quantum point contact in the FQH regime like in Fig. ^was realized by 
Saminadayer and Glattli ^H]- Near the potential V{x,y) of the constriction the filling factor 
was u = ^. For small currents Ib (weak backscattering limit) the measured shot-noise jl5M16j 
in the current / let them infer a fractional charge of |e for the particles that are tunneling. 

The regime of strong backscattering Ib ^ I was investigated by Griffith et. al. ^U] ^-nd 
showed behavior as expected: If the potential of the barrier is sufficiently high the system is 
effectively divided into two halves and only electrons can tunnel between them. 

Just recently Chung et. al. |2Uj made experiments at lower temperature on a system of two 
quantum point contacts. The quasiparticles flowing along the edge were reduced in density 
by transmitting them through the first quantum point contact prior to hitting the second 
one. Thus being in a regime where quasiparticles arrived "one by one" at the constriction, 
single quasiparticle tunneling could be investigated in the absence of correlations between the 
particles. The measurements showed that single-quasiparticle tunneling was only observed 
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Figure 1: Proposed geometry for quasiparticle tunneUng between one-dimensional edges in the 
FQH regime. The source-drain voltage Usd causes a current / through the sample. 
The edges are bent due to an external potential V{x, y). Tunneling of quasiparticles 
is expected to occur where the edges are close to each other. The tunneling current 
Ib causes fluctuations in the transmitted current / which is measured. 
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if the temperature was sufficiently high (73 mK) while at lower temperature (23 mK) only 
electrons were found to tunnel even at quite transparent constrictions. This very surprising 
result is supported by recent calculations by Kane and Fisher ^\ founded on the chiral 
Luttinger liquid theory. Their calculations show "strong evidence" for only electrons being 
transmitted at T = and their proposed explanation is an Andreev reflexion. 
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3 Theoretical preparations 

In this section the preparations for the numerical calculations on systems in the fractional 
quantum Hall regime are collected and derived. In section IXTl the system to be treated will 
be introduced and the necessity of exact diagonalization will be pointed out. The subsequent 
sections deduce the single particle basis used for the numerical calculations. In this context 
the periodic boundary conditions will be explained and a physical interpretation for the 
appearing phase factors will be given. In section I.S.2I using these boundary conditions a 
possible realization of an in-plane electric field will be discussed. 

3.1 The model of the fractional quantum Hall system in this work 

The model of the fractional quantum Hall systems treated in this work make the simplification 
of electrons to be confined in a two dimensional plane. Perpendicular to this x-y-plane 
there is a homogeneous magnetic field B = Bcz- The electrons interact via a repulsive 
interaction potential Vint, for which we will use the Coulomb interaction as well as a short 
ranged interaction and in some cases they are subjected to an external potential Vext which is 
used to model a constriction inside the system. Additionally, there can be an in-plane electric 
field used as a driving force for a current. 

Accordingly, the many particle Hamiltonian has the following form 

^ = ^ E + E ^-*(^"*^) + E -^^y (1) 

i=l i=l i<j 

Here, H is the kinetic momentum that incorporates the vector potential caused by the mag- 
netic field and the in-plane electric field. We will calculate in the limit of high magnetic field, 
B ^ oo. So the spin of the electrons is assumed to be polarized and there is no spin degree 
of freedom in the Hamiltonian. The states of electrons in a magnetic field are known to be 
quantized in macroscopically degenerate Landau levels En of energy £"„ = TiuJc{\ + n) with 

I , — eB 
'- m 

Here we are interested in the case of a partly filled first Landau level, mainly in the filling 
factor V = = where Ng denotes the number of states in the lowest Landau level. In 
the limit of high magnetic fields the energy is dominated by the energy of the Landau level 
quantization. The degeneracy of the Landau levels forces us to use exact diagonalization to 
account for the interaction and for external potentials. Thus it is reasonable to determine 
the eigenbasis of the single particle kinetic Hamiltonian H^in and treat the interactions and 
external potentials by diagonalization in the degenerate space of the lowest Landau level. 

This basis will be derived in the following sections. 

3.1.1 One electron in a magnetic field: Harmonic oscillator 

Consider electrons that are constrained to move in the x-y-plane and which are exposed to 
a homogeneous magnetic field B = Bcz parallel to the z-axis. The vector potential can be 
chosen in Landau gauge as 



A{r) = B X ell, 



(2) 
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where r = {x,y) is the coordinate in the plane. 

The Hamiltonian describing one electron in this system is 



d 

with Hx = —ifi'TT- 
ox 



d 

Hy = —ih— — h eB X 
dy 



The operators Tlx = Px + eAx and Hy = Py + eAy are the kinetic momenta of the electron, P 
denotes the canonic momentum. A natural energy scale of the system is given by the cyclotron 

frequency = the typical length unit is the classical cyclotron radius Iq = ^J^- By 
rescaling the kinetic momenta to these natural units two operators Q and S can be defined 
as 

Q = -r^^y (4) 
S = , ^ Ux. 



Their commutator evaluates to 

[Q,S] = (5) 

which is the canonical commutation relation. 

Following |21j . it is possible to express the Hamiltonian (jH)) using Q and S. Knowing 
the commutator Q one recognizes this to be the problem of an one-dimensional harmonic 
oscillator, whose spectrum consists of equidistant discrete energy levels, commonly known as 
Landau- levels. 

H = huJciQ^ + S^) (6) 
En = huJci^+n); nGN° 

In an analogous manner to the approach for solving the usual harmonic oscillator, raising and 
lowering operators respectively a can be defined as 

at = -^{Q-iS) (7) 



With aid of the commutators [Qi'S'^J = 2iS and [5,(5^] = —2iQ following from © the 
action of or a on an energy eigenstate of the system can be verified to raise and lower its 
energy by an energy quantum huJc- 
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3.1.2 Magnetic translations 

Once chosen the Landau gauge, the Hamiltonian @ given above commutes with the momen- 
tum operator Py in y-direction, since 11^ and liy do. So it is possible to find a complete set of 
common eigenvectors \n,py) to both operators H and Py, where the label n denotes the energy 
eigenvalue and py the eigenvalue of the momentum. So we have H\n,py) = hu!c{n + ^)\E,py) 
and Py\n,py) = py\n,py). Apart from the raising and lowering operators for the energy there 
must be a momentum-shift operator B^{dp) affecting the eigenvalue of Py like Py ^ Py + Sp. 
This implies the following commutator of Py and 

PyB^{6p)\E,Py) = (Py+Sp) B\5p)\E,Py) (S) 
B\6p) Py\E,Py) = B^{6p)Py\E,Py) 

^[Py,B^{6p)] = 5pB\5p) 

A simple momentum translation operator of the form B\6p) = exp(^(5p y) does however not 
commute with the Hamiltonian, which is due to the non vanishing commutator between y 
and lix- But a generating operator = y — i-^dx commutes with Yix as well as with Hy 
(compare Equ. (jSJ) and is thus compatible to H. A finite transformation constructed from 
this generator is a momentum shift by 5p accompanied by a coordinate shift in x-direction 
by X, = 'fSp. 

6t = (9) 



B\5p) = exp(— (5p 6^) 



exp [^5p y^ exp {j^dx 



tx{Xp) := exp ( "2 Xp y ) exp ( Xpdx 





Due to the coupling of x coordinate to the y-momentum, we define tx{^^) := B^{5p) as the 
"magnetic translation operator in x-direction" . Since tx has the demanded properties of (jH)) 
it is the (continuous) raising operator for the eigenvalue of Py. This continuous symmetry of 
the system causes an infinite degeneracy of every Landau level. 

Apart from the magnetic translation tx{Xp) in x-direction, the operator Py can be used as 
a generating operator of an ordinary translation in y-direction, that as well commutes with 
H 

ty{Y) = eMYdy). (10) 



3.1.3 Basis in direct space 

By virtue of the operators a, ,t{6p) the eigenbasis can be determined in direct space. Initially 
the ground state can be calculated by applying the lowering operator. We are looking for the 
simultaneous eigenvector of H and Py for the eigenvalue in both cases. This eigenvector 
has to satisfy 

a|0,0) = (11) 
Py\0,0) = ^ 0o,o(2;,y) = 0o,o(a;). 
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The second equation tells us that (pQ^ is only a function of x. Now using the direct space 
representation of the lowering operator, we gain a differential equation for (/)o,o(^) 

(nj,-m^)|o,o) = (12) 

{—ihdy + eBx + hdx)4>{x) = 

1 

^(t>oflix,y) = —==exp{-—^), 



which, as a function of x, is normalized to 1. By applying B'{py) = tx{-^Py) one obtains 
solutions with a different momentum py in the same Landau-level. This leads to the explicit 
form of eigenfunctions in the lowest Landau-level 

\0,Py) = U^py)\0,0) (13) 



1 fi \ f-{x + Xp) 



with X„ = -^py 
n 

which are normalized to a delta distribution (like plane waves are). The eigenfunctions of 
higher landau-levels can be obtained by the repeated application of the rising operator 
which results in an additional factor consisting of a Hermite polynomial Hn{^^j^) of degree 
n where n denotes the Landau-level. 



3.1.4 Generalized periodic boundary conditions on the torus 

The aim to describe an infinitely big system in a finite basis that can be handled for numerical 
diagonalization, can be achieved only if the system complies with certain symmetries, which 
reduce the degrees of freedom to a finite number and thus allows for a description in a finite 
basis. Here, the system will be required to have a discrete translational symmetry where the 
periodicity is given by the size of a cell {x, y) S [0, Li] x [0, L2]. Thus all physical observables 
have to return to their original value when proceeding to another cell. 

To begin with, the number of electrons in the x-y-plane must be reduced to a finite number 
in order to have a finite number of degrees of freedom. The idea is the following. Only a 
small number of electrons is put into one unit cell [0, Li] x [0, L2] and treated independently. 
All neighboring cells only contain images of these electrons at exactly the same position 
with respect to this cell's boundaries. The intention is to model a system that spreads out 
infinitely and doesn't need any confining potentials which would cause boundary effects. 
Since the interaction of electrons in distant cells is somewhat smaller than that between 
electrons in the same and directly neighboring cells, the mirroring of electrons is assumed 
not to destroy the desired effects. However, to gain results of a real infinite system with 
infinite degrees of freedom, a careful analysis for the size dependence of the observed effects 
would have to be made. A sketch of the repetition of the unit cell is found in figure |2j The 
equipotential lines are intended to depict a potential which is periodically repeated in every 
cell. Putting this more formally, we have a finite number of electrons in the system. Due to 
the cell wise periodicity of all potentials in the system, we know by Bloch's theorem, that 
the eigenfunctions of the one particle Hamiltonian will be simultaneous eigenfunctions to 
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Figure 2: Periodic tiling: The physical content of the unit cell -i.e. external potentials 
(equipotential lines in the figure)- is repeated over the entire system. A finite 
number of electrons (blue circles) in this cell is treated independently, the electrons 
in neighboring cells are mirrored. Although only depicted for the y-direction, the 
periodicity holds in x-direction as well. 
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the translations that do a shift by one cell size. Since here we have a non-periodic vector 
potential in the kinetic part of the Hamiltonian, instead of ordinary translations we have 
to use the magnetic translations tx{Li) and ty[L2) from section [3.1.21 in order to commute 
with iiffcj^.In complete analogy to Bloch's theorem, these translations are unitary operators 
which have eigenvalues of modulus 1, which is a phase factor the wavefunction picks up 
upon application of tx{Li) or ty{L2) respectively. Thus, we arrive at generalized (because of 
magnetic translations) periodic boundary conditions 

tx{Li)4){x, y) = ex.-p{ia)(t){x, y) (14) 
ty{L2)cl){x,y) = exp{i^)4>{x,y), 

where a and (3 fix the eigenvalue for the respective translation. 

The above procedure is of course only possible if additionally and ty commute. Generally 
this is not the case, since their commutator evaluates to 

[^y + dx,dy] = -^GC (15) 



ty{L2)tx{Li) = tx{Li)ty{L2)exp (^■^LiL2y 



The appearing phase factor exp (^p-LiL2 j is the phase factor known from the Aharanov-Bohm 

effect 1^ . The connection between these to phase factors can be understood in the following 
manner. One has to keep in mind that the magnetic translations perform a coordinate shift 
with the additional property of transforming eigenstates of -fffcm to eigenstates with the same 
energy. In the chosen Landau gauge, the vector potential A(x, y) is linear in x and independent 
of y and thus ^(0, 0) = 0. Assume we had an eigenvector of i^fcm ^-^d want to construct 
another eigenvector which is shifted by {a,b). In a first step we can just shift the reference 
frame by {—a, — b) which changes x ^ x + a, y— + 6 

Hkin\x^x+a,y-^y+b "^oi^ + a,y + b) = Eo'^oix + a,y + b) 

^ V ' 

H' 

"kin 

The coordinate shift in Hkin only affects the vector potential which accordingly changes by 
A ^ A' = A + eBaCx- The new Hamiltonian can be written in a more complicated form 

exp(- / AdT)H'^,^ = Hkinexpi- A dr) (16) 

" Jc{a,b) Jc{a,b) 
%€■ f — » 76 f — » 

Hkinexp{— / Adf)'^o{x + a,y + b) = Eexp{— / A df)'^o{x + a,y + b) 

Jc{a,b) Jc{a,b) 



The phase factor compensates for the offset eBa in the vector potential and it is the same 
factor that appears in the explanation for the Aharanov-Bohm effect. Here c(a, b) is a contour 
that connects the point (0,0) with (x + a,y + b). Since the vector potential is not rotation 
free, ^a,b{x^y) is a functional of c and depends on the chosen path to reach the endpoint. 
Therefore, the exponential factor in El also depends on this path. 
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Figure 3: Aharanov-Bohm effect. First applying tx{Li) and then ty{L2) shifts an eigenvector 
along the blue path, the reverse order moves it along red path. The solutions in 
the point (Li,L2) differ by the Aharanov-Bohm phase related to the magnetic flux 
enclosed by the two paths. Commuting operators and ty imply an integer number 
of enclosed flux quanta (green circles). 



This procedure can be regarded as an alternative derivation of the magnetic translations. 
Comparing Equ. Eland Equ. ©, we obtain tx{Li) from the first equation by choosing c the 
path that is parallel to the x-axis from (0, 0) to (x + Li, y). 

Analogously we can construct the translation ty, where the phase factor is trivial. 

By successive application of tx{Li) and then ty{L2) or first ty{L2) and then tx{Li) to 
an eigenvector, move the solution around the unit cell from point (0,0) to (Li,L2) on two 
different traces. These two traces are sketched in figure 01 The Aharanov-Bohm phase by 
which the two solutions differ is just the phase factor in Equ. I15L This phase factor can 
be written like in (|17() and -by Stokes' theorem- be related to the amount of magnetic flux 
enclosed by the contour. 



To return to the commutation relation of tx{Li) and ty(L2), commuting operators demand 
a phase factor which is a multiple of 27r in Equ. 1151 Thus the wavefunction must pick up a 
phase factor 27riVs with Ng integer. Writing the wavefunction as a function of the complex 
variable z = x + iy, this corresponds to a complex function with A'^, zeroes in the unit cell. 
This number Ng on the other hand is the number of flux quanta $o = f piercing the unit 
cell, which can be seen from the expression of the total flux <I> through the area L1L2 of the 
unit cell 

$ = L^L^B = Ng- = Ns^o (17) 

e 

L1L2 = 2ttNsII, 

where the definition of the magnetic length Iq = was used. 

From the second equation in (|14jl . keeping in mind the form of the wavefunction in (|13j) . 
the quantization of the y-momentum follows 

Pj = ^-rJ + -rP J e z. (18) 

1-12 -1^2 
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Keeping in mind that the magnetic translation tx{Li) = B'^{pLi) = B"^ (Ngfi^) has the 
meaning of a raising operator for the eigenvalue py and provided the periodic boundary 
condition ()14() . this is an identification of the states \n,py) and \n,py + kNgfi^), where k is 
an integer. Thus we are looking for a linear combination of those identified vectors satisfying 
the periodic boundary conditions. These will be named \n,j)prb and can be constructed from 
the non-periodic solutions of (|13|1 by summing up all identified wavefunctions 




To obtain the wavefunction in direct space, and tx have to be expressed in this represen- 
tation, which then results in 

, , , 1 ,Xi + kLi , . / (x + Xj + kLi)"^ \ 
Mx,y) = 2^exp(z y-z/ca)exp( 1(20) 



fceZ 

X H, 



X -\- Xj + kLi 

Here Hn are the Hermite polynomials. This wavefunction is normalized to unity upon inte- 
gration over the domain of the unit cell. 



3.1.5 Physical meaning of a and /? 

The two parameters a and (3 were introduced in section 13.1.41 to fix the eigenvalues of the 
magnetic translation operators according to Equ. p4|l. A physical meaning of magnetic 
fluxes that pierce the torus can be attributed to the parameters a and (3 as shown by Tao 
and Haldane [23]. The parameters introduced in their paper can be shown to be equivalent 
to ours by the following consideration. 

The periodicity in the tilings of one unit cell can alternatively be regarded as an identifi- 
cation of the left border of the unit cell with its right one and the upper with the lower one. 
This imposes a torus topology. Because the x- and y- axis differ only by the selected gauge, 
it is sufficient to focus on one parameter, here (3. If we apply a unitary transformation 

U{y) = eM-ij-y) (21) 

on both, the Hamiltonian and the periodic eigenvectors, the new eigenvectors \k) satisfy 
"simple" boundary conditions with respect to ty{L2), like 



ty{L2)U{y)\k)=U(^ 
\k) 



(22) 
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Figure 4: The parameter (3 can be interpreted as the magnetic flux of a solenoid parallel 
to the y-axis, if x = and x = Li is identified by the periodic boundary condition. 



Therefore, both borders can be identified by means of ty{L2) like (x, y) ~ (x, y + L2). The 
transformed Hamiltonian gains an additional term residing in the momentum operator liy, 
which transforms according to 

iiy = U{y)YlyU\y) = Iiy + ^p. (23) 

The additional term proportional to (3 can be seen as belonging to a vector potential Ap = 
^Sy. This again is related to a flux of = |/3 passing through the torus as sketched in 

figure 1^ If now (3 is varied with time, an electric field Ey = —Abeta = ~^r^^y is induced in 
the plane along y-direction. 

To stress, that this interpretation is reasonable, the influence of (3 on the basis states from 
Equ. ()2U() (which still have to be multiplied by U{y)) is reviewed. An electron in the state 
\k) is localized around —Xk^p in x-direction. Since Xk^p depends on (3, the electron moves in 
x-direction is f3 changes adiabatically. The velocity of this motion is 

(24) 



Ey 

B 

which coincides with the classical E x B drift. 



L2L1 

2ttN, h ^ 



'0 



3.2 FQH system with an applied electric field 

The parameters a and /?, formally introduced in the periodic boundary conditions in section 
I3.1.4( turned out to have a physical meaning of magnetic fluxes as shown in the previous 
section. It was shown that a time dependent change of /? can be used to generate an electric 
in-plane field in the system. The same is true for a since x- and y-direction only differ by the 
selected gauge. In this section we will assume to describe a system in the basis of states in 
the lowest Landau level (n = 0) from equation (|19j) and use /? for creating an applied electric 
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field in y-direction. In time dependent perturbation theory the time evolution of the system 
will be calculated. 

A homogeneous magnetic field B = Biz and an electric field e = ecy in both can be 
represented by a vector potential depending on space coordinates and time 

r)A 

B = V X A 
= 
A{r,t) = {Bx-et)ey. 



The part of A describing the magnetic field was chosen in Landau-gauge which makes it 
possible to use periodic boundary conditions (|14)) with fixed values of a and (3, say 0. The 
electric field in y-direction can be achieved by a homogeneous term of A linear in time. 
Therefore the many-particle Hamiltonian becomes time dependent and obtains the following 
form 

Ne . Ne 

= T.^(^lx + ^JAtf)+T.^^nti\rJ-n\) + Y,ye.tirj) (26) 
i=i j=i 
d 



dxj 



d 

Ilj^y{t) = —ih-^— + eB Xj — ee t. 

Here, Hx and IIj^ are the kinetic momentum operators, Vi,if(|r|) is the interaction potential 
and Vexti^^ is an external potential compatible with the periodic boundary conditions, needed 
to model a constriction. 

It is possible to isolate the Hamiltonian's time dependence from (|26|) in a time dependent 
parameter (3{t). The second parameter a influences the momentum in x-direction. Both 
parameters — a and /3 — can as well be cast into modified periodic boundary conditions 
by a gauge transformation Ua^p{x,y) = exp(— ■|(xa + y(3)). This will yield the generalized 
periodic boundary conditions from section ()3.1.4() . Either the boundary conditions are kept 
fixed and the parameters are regarded as parts of the Hamiltonian or the Hamiltonian will 
have a = /3 = and the periodic boundary conditions will catch up a phase factor of exp(ia) 
upon a magnetic translation one unit cell to the right or a factor of exp(i/3) when going one 
cell up, respectively. However, now we will keep the boundary conditions fixed and concern 
the parameters as parts of the Hamiltonian 

Ne Ne 

= Y.^^^lj,- + ^h,y) + Y.^-oui{P,-rk) + Y.Vext{f,) (27) 

j=l j^k j=l 



d G fx 

^p,j,y = -i^i^ \- -B Xj + —13. 

dyj c L2 
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Here, Li and L2 again denote the extent of the unit-cell. The time-dependent problem to be 
solved thus is the solution of the Schrodinger-equation of the time-dependent Hamiltonian 

H(') = Ha=o,,=-^f (28) 

3.2.1 Time evolution in time dependent perturbation theory 

The time dependent perturbation theory allows to calculate the state's time evolution in the 
case where the time dependent terms in the Hamiltonian can be considered to be small with 
respect to the unperturbed Hamiltonian Hq, where 'small' compares the typical energy of the 
perturbation to the spacing of the eigenenergies of Hq. In the adiabatic limit the perturbing 
term is switched on smoothly in an infinitely long period of time while the total change of H(t) 
is yet finite. Kato has shown that even in the case of crossing levels the time evolution 
will follow the stationary eigenstates of H{t), for fixed t, except for a phase factor (see next 
section). Applied to our system, when starting in an eigenstate of H{0) = i?Q,=o, /3=o in t = 
the system will follow this state on changing (3 adiabatically thus staying in the corresponding 
eigenstate of the operator H{t) = H^^q^ ^(j) for all t. Now we want to consider in addition 
the phase factor which is lost when using the adiabatic theorem. This phase factor is crucial 
for the task of constructing the time evolution operator of the system from its eigenstates. 

For every time t (and consequently every value of a and /3) we will therefore consider a 
spectral decomposition of H according to 

Ha,l3 = ^ \ka,l3){ka,f3\^k,a,f3- (29) 
k 

Starting in an eigenstate |A;o,o) for t = 0, the adiabatic evolution would correspond to following 
the state |fco^/3(()) for all t. If however the electric field is not switched on adiabatically, we can 
still use this time dependence as a starting point to govern the true time evolution by means 
of perturbation theory using the adiabatic states as a basis. Since the basis is complete (in 
the subspace of the lowest Landau level), an arbitrary state in the lowest Landau level can 
be written as 

^{t) = 2_,afc(t)exp(-- / Sfc,a,/3(t') ^^^0 I^a,/3W) (30) 
k ''^ 

Ofc(t)exp(- / Eka^p, dp) \ka,(s). 

k ^2ee Jo 

Using this ansatz in the time dependent Schrodinger-equation we obtain the following result 
(we omit here the dependence on a for the sake of simplicity) 



ih^-^ = Hmit) (31) 



^(«fc(i) + ^afc(t)£;fe^^(i))|A;^(t)) -Fafc(t)^|/c^(())^ x 

X exp(-- / £^fc,/3(t') dt') 
Emit) = ^afe(t)£;fc,^(t)exp(-- / ^fc,/3(t') c?i')l%W>- 

I, -^0 
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By projecting on the state (fc/3(t) I and using the basis' orthonormality we arrive at a differential 
equation for the coefficients afc(t) 

ihdk{t) = i?i^a;(t)(A;^(i)|^i^^exp(-^ / (^/,/3(t') - -^fc,/3(t')) ^*')- (32) 

In the first place, the sum over / includes off-diagonal elements as well as diagonal terms. But 
it turns out however that only the diagonal terms contribute, which can be shown with help 
of Hellmann-Feynman's theorem [201, by calculating (A^/jj^^p^ = ""^^(^l^f^- 

Hellmann-Feynman's theorem is applicable on parameter dependent eigenvectors of a Her- 
mitian operator dependent on the same parameter. Normalized eigenvectors assumed, we 
have -^{kpllfs) = 0, from which follows —^jj^llfs) = ~(^/3l^^- case we yield for Hf^ 

{kp\Hp\lp) = Ek,p6k,i (33) 

Here S^^i denotes the Kronecker symbol. With H27() the derivative of can be expressed as 

~ ^P'j^y restrict the state to lie in the sub-space of the lowest Landau- 

level, the matrix elements of the kinetic momentum operators 11 are zero. This is because 11 
can be expressed as a linear combination of and a. It will be shown explicitly in section 

Em 

From equation (|33|) for k = I we would conclude that the eigenenergies should be constant 
with P and thus with time. This however is only true for homogeneous systems. The reason 
for that is the projection to the lowest Landau-level. In the language of perturbation theory, 
the eigenstates are given to lowest order, while the eigenenergies are calculated to one order 
above. A more precise treatment will be given in section EiS For the moment it is enough to 
assume E/^^p to depend on (3. 

So we have as a final form the differential equation for the coefficients ak{t) 

iMkit) = ak{t)ih{kp(^t)\^^^^ (34) 



akit) = afe(0)exp(^y" (^>(t')|-^^f^*'^ 

This is a phase factor as the argument in the exponential function in 1)34(1 is an imaginary 
number, since {k\^ = -^\k) = -{k\^ . 

Given an arbitrary initial state of the system, we can calculate the time evolution of the 
system by projecting the state onto the eigenstates for t = which yields the coefficients 
afc(O). These then evolve according to ()34|) . The state of the system for the time t is the 
superposition of these contributions given by equation (jSOJ. Comparing this result to the 
adiabatic limit in the next section, here we have an additional phase factor of equation (|34|) . 

3.2.2 Adiabatic limit of time evolution 



To appreciate the results from the previous section, they will be compared to the adiabatic 
limit of the time evolution. In the limit of an infinitesimally weak electric field the time 
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evolution of a state can be calculated by means of the adiabatic theorem, first used by Fock 
and Born, later generalized by Kato ^7\. Here we use it in a similar way as in |2H]- This 
theorem facilitates the calculation of the time evolution of a system described by a parametric 
dependent Hamiltonian, where the parameter is varied infinitely slowly by a finite amount 
or — put differently — the total change of the Hamiltonian from t = — oo to t = is finite. 
Following the nomenclature in [2ZI the time dependent Schrodinger equation is rewritten by 
rescaling the time variable hy t = ts. This makes it possible to separate the parameter s 
running through a finite range from r which goes to infinity 

indt^{t) = H{t)-^{t) (35) 
^itids'HAs) = TH{s)-^r{s). 

The aim is to arrive at an expression for the time evolution of the system in the limit of 
T oo, namely ^{s). The adiabatic theorem as proven in the paper of Kato states that for 
an eigenfunction 0(0) of H[{)) there exists the following approximate expression for the real 
time evolution formally described by the unitary operator Vr{s) 

Vr{s)(t>{<d) = exp E{s) ds^ (j){s) + 0{t-^) (36) 

where if(s)(/)(s) = E{s)(p{s) 

In the limit of r — > oo this expression becomes exact, indicated by the vanishing deviation- 
term. 

Applied to the problem of an additional electric field, we can determine the time evolution 
for infinitely weak electric fields. According to equation (|28|) r can be identified with the 
parameters introduced earlier as 



P = -^t. (37) 




Using (|57j) in equation we arrive at 

^{(3) = lim exp(-^ 1'^ E{f3)dp) m- (38) 

e^O \L2ee Jq J 



So the phase factor in equation H34() would be lost by this adiabatic treatment. 
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4 The homogeneous system 

4.1 Laughlin's wavefunction, interaction and correlations 

In this section some well known features of the homogeneous system are collected and partly 
adapted to our geometry, since they are crucial for understanding and motivating the quasihole 
excitations. 

For the circular gauge of the vector potential in the absence of any background potential, 
Laughlin has proposed in j3| a Jastrow-type variational wavefunction for the case of filling 
factors 1^ = ^ where m is an odd integer. 



Picking one electron's coordinate, say zi, while keeping the others fixed, the wavefunction 
exhibits (iVg — l)m zeros, m of which appear whenever zi approaches another electron's 
position. This m-fold zero produces strong correlation holes around each electron which in 
turn reduces the Coulomb energy. Since the number of zeros of the wavefunction is fixed by 
the number of flux quanta through the system, the only freedom is where to put these zeros. 
In equation H39() the maximum number (Nf. — l)m of zeros is used to create the most effective 
correlation holes. These are the reason for this ground state to be separated by a gap from 
the system's excited states. Numerical comparisons showed for different repulsive interaction 
potentials that this trial wavefunction has a big overlap with wavefunctions found by exact 
diagonalization [31 Ej. Independently Yoshioka, Halperin and Lee [SJ found these correlations 
in their numerical work. 

Along with the ground state Laughlin also found the elementary excitations of the system 
to be quasielectrons and quasiholes. Laughlin mapped this system by an analogy to an one 
component plasma and thus identified these quasiparticles to have fractional charge it— e. 
These excitations are created whenever the filling factor u deviates from its value — : If it 
increases, quasielectrons are created; on decrease there are additional quasiholes. The finite 
amount of energy e± needed to create one of these particles makes the system incompressible, 
because an infinitesimal change in the system's area (which of course changes ly) causes 
quasiparticles to be created each of which rising the energy by e±. Thus the compressibility 
is infinite. This incompressibility will be verified in section l4.1. 51 

4.1.1 Short-range interaction 

Although the overlap of Laughlin's wavefunction with the exact one is high for Coulomb 
interaction it is not an eigenstate of the system. Haldane and Rezayi |34j treated the interac- 
tion by means of introducing pseudopotentials with different ranges. They showed that the 
ground state mainly depends on the pseudopotential parameter of shortest range and that it 
is quite robust against changes of the other parameters. They found out that the Laughlin 
wavefunction is an eigenfunction in the case of a short range interaction which only has the 
one nonzero pseudopotential parameter for the shortest distance. This treatment was specifi- 
cally used on the spherical geometry and relied on conservation of angular momentum, which 
is not true for our system. 

This short range interaction was generalized by Trugman and Kievelson j^H] who wrote 
down an analytical form for it without resorting to angular momentum conservation. Al- 




(39) 
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though they used the open plane geometry we can follow their idea here and show that it 
is also applicable to rectangular geometry. In the paper cited above, an expansion of the 
interaction potential is considered as Vint in terms of its range b as 

Vnt = J^c,-62JV2^,5(f). (40) 

j=0 

Going to the Fourier-space, this expansion is seen to be a Taylor series of a symmetric function 
in \q\ (only even powers of \q\ appear), which is true for every real valued isotropic potential 
which has an analytic Fourier transform in g = 0. The requirement for the potential to have 
no singularity in g = is fulfilled if the potential has a finite mean value, since the Fourier 
component V{0) = dr rV{r) is the mean value of the potential. Thus, if V{q) is analytic 
in g = 0, it is expandable into a power series in |gp and we can determine the coefficients 
Cj of Equ. H4U|) . This for example is true for a Yukawa-potential. The Coulomb potential 
however does not have this property. 

The first term co(5(r) in ()4U() vanishes for spin-polarized antisymmetric wavefunctions. The 
leading term in the limit of small ranges b thus is 

Vshortif) = b'vHir). (41) 

This interaction was shown [221 reveal Laughlin's wavefunction as the unique ground state 
for V = ^. Therefore, using this interaction instead of the Coulomb potential, we can assume 
the Laughlin-like correlations in the wavefunction, which are manifest in the relative zeros of 
two electrons' coordinates, to be more pronounced than for Coulomb interaction. 
The Fourier transform of this short-range interaction is 

Vshort{q) oc -|gf . (42) 

In the case of inhomogeneous systems (section ^ the work of Krause-Kyora [201 showed 
that using this short-ranged interaction is advantageous due to the absence of the long-range 
part of the Coulomb potential which causes oscillations in the density profile even far off the 
actual constriction potential. Therefore, we will make use of this interaction in parts of the 
following calculations. The term "short-range interaction" and "hard-core interaction" will 
be used equivalently in what follows. 

Apart from introducing Laughlin's wavefunction in the next section, the formal reason for 
this interaction to make Laughlin's trial wavefunction an eigenstate will also be investigated. 

4.1.2 Laughlin's wavefunction in rectangular geometry 

Haldane and Rezayi [22 ported Laughlin's wavefunction the rectangular geometry with peri- 
odic boundary conditions and they investigated impurity effects on this ground state. From 
this paper we will cite the many particle wavefunction here since it will be used in section 
14. 1.31 to show that the expectation value of the short-range interaction vanishes for this state 
and also to construct the quasihole creation operator in section 14.2.11 

The starting point is a Jastrow-ansatz which implies a many-particle wavefunction that 
separates into a product of a relative- and center-of-mass-function, where the relative part 
factorizes to a product of functions /(zj — Zj) each depending only on the difference of two 
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particles' coordinates zi^Zj with zi = Xj + iyi. The center of mass coordinate is given by 
Z = Zj. Thus we can write 



^>{z^,...,ZNe) = J[f{z,-Zk)FNAZ)eiCp\-TTj-f-Y.''^- (^3) 

j<k 




The wavefunction F/v^ (Z) exp(— tt ^ £ Xli^?) ^^^-s the same form as a single-particle wave- 



function of an electron in a magnetic field containing m = ^ flux quanta. Its explicit form 



122] is 

FN^iZ) = expiKZnfl^i^TT^^^lt^^ (44) 

where exp(— zA'L2) = (—1)^° exp(i/3) 
and exp ^ Z*^ = (-1)^= exp(ia + KLi), 

where i?i is the odd elliptic theta function of first kind |2H|- It is important to note that 
"diiz) oc z for small \z\ such that F^^ has m zeros. The solution of (|44() is uniquely defined 
by giving the real wavevector K and these m zeros Zj/. The number of linear independent 
solutions of (|44() were shown [22j to be equal to the number of zeros Zi,, thus there are m. They 
cause an m— fold degeneracy of the ground state of a homogeneous system. Alternatively, 
according to Tao and Haldane [21], this degeneracy can be regarded as originating from the 
translational invariance of the center of mass wavefunction. 

The function f{z) in|121has to be chosen such that the periodic boundary conditions are 
fulfilled and they have to be odd functions in order to obey the Pauli principle. Like in the 
open plane, demanding an m-fold zero whenever two electrons approach each other leads to 
the single solution 



„ / -i(z*-zl) L^\m t^NsTJ^I.x'] 
^Gs{zu...,ZMe) = n^i(^ r l^r) ^^»(^)^^P( LL ^^^^ 

j<k 2 2 12 



Fixing all but one electron's coordinate and counting the zeros of the wavefunction as a 
function of the last free coordinate, we clearly should find zeros within the unit cell in 
accordance with the Aharanov-Bohm effect, as stated in section l!j.l.41 It is important to note, 
that we will find an m-fold zero whenever approaching one of the fixed electrons' coordinate, 
giving (A^e — l)m zeros in total. This is seen from ()45() . The remaining m zeros can be found 
in the center of mass wavefunction (see Equ. (|44() ). Their positions obviously depend on the 
positions of the other electrons. This is an important point to make, because it shows that 
we cannot fix these zeros with respect to one of the electrons at a certain point. Instead, if 
we want to fix a zero of the wavefunction at some given point with respect to one electron's 
coordinate independently of the other electrons' positions, we have to insert an additional 
flux quantum to gain the freedom of localizing the respective zero arbitrarily. In section 14.2.21 
this will be used to create a quasihole. 
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4.1.3 Vanishing short-range interaction for Laughlin's wavefunction in rectangular 
geometry 

In the previous section the equivalent to Laughhn's wavefunction in a system with periodic 
boundary conditions was cited. What Trugman and Kievelson [2^1 did for the open plane 
geometry is possible to show for periodic boundary conditions as well: This wavefunction 
yields a vanishing expectation value {Vghort) for the short-range interaction, Vghort given by 
Equ. (|4H) . Furthermore, also the trial wavefunction for a quasihole which will be given in 
Equ. (|59() minimizes this interaction to zero. 

To proof these two statements, we can consider a trial- wavefunction of the following form 



^{zI,...,ZnJ = Yl^{Zi- Zjrf{zi,...,ZNj 



(46) 



i<j 



Here 'd{z) = 'di{—iirj^\ij^) is used as a short hand and f{zi, . . . , zn^) is a symmetric function. 
7?i is the odd elliptic theta function of first kind. In the case of Laughlin's ground state 
wavefunction the function / is given by / = F/v^ (Z) exp(— tt j^^£^ X^i^D' with Fn^{Z) 
being the center of mass function given in (|44|) . If we take the trial wavefunction of a quasihole 
excited state, / = Il^^-^'d{zk — zq)F]^^j^i{Z + ^) exp(— tt^^^ Ylii ^f) {^^ can be seen from Equ. 
(|59|) ). The exact form of / however does not matter here; only the fact that it is 2 times 
differentiable and that it doesn't have any singularities will be used in the following. 

Writing the coordinates z = x + iy as r = xe^ + ye^, the expectation value of the two body 
interaction can be calculated as 



Vshort) 



/ dVi...(iV^J^'(zi,...,z^e)|2 ^y(ri-rj) 

^^e(iVe-l) f ^2^^... ^2^^^ | *(.l , . . . , .^e) P (H - r^) . 



(47) 



Introducing new coordinates r_ = ri — r2 and r_|_ = ri + r2, this integral can be recast into 



{Vshort) 



I 



X / d^r. 



d^rs . . . d^VN, / r+ x 



f(zi,...,ZNe)Y\_^i^i - ^jY 

i<j 



(48) 



where the short-range interaction potential was substituted. The last integral I can be inte- 
grated by parts twice, resulting in 

2\ 



I 



/(Zl, . . . ,ZAre) JJl?(^ 



(49) 



f{zi,...,ZNe)Y[^{ 



Zi Zj 



i<3 



r_=0 



The boundary terms vanish for both steps of partial integration. The boundary term after 
the first integration contains a factor V5(r_) with r^^x = or r^^y = ibL2 as limits. 
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After the second integration, the factor (5(r_) in the same Umits makes the boundary term 
vanish. Also the last expression ()49() vanishes, because after differentiating we can factor out 
\'&{zi — 2:2)1^"^"^ that behaves like \zi — Z2p'"~^ if r_ — > 0. Thus, the positive semi-definite 
interaction energy is zero for both — the ground state and the quasihole state — and, since this 
is the only contributing operator in the Hamiltonian (apart from the frozen kinetic energy), 
those states must be ground states of the homogeneous system. 

Going to higher powers in j in Equ. H4U|) and applying the same argument as before, the 
expectation value vanishes as long as j < m. However, since we are interested in a filling 
factor V = ^ = the order j = 1 in the short-range interaction is already sufficient to 
force a threefold zero to reside on every pair of relative coordinates in the ground state. 
The reason for that is the Pauli principle which demands an odd number of zeros on every 
relative coordinate. So, the interaction potential given in (|4H) has the property of producing 
Laughlin's wavefunction H45|l in the case of a filling factor 1/ = |. 

4.1.4 Projection to the lowest Landau level 

So far we have seen trial wavefunctions for the homogeneous system. Now the many particle 
Hamiltonian given in Equ. ^ of section 13.11 shall be treated by numerical diagonalization 
first in absence of external potentials Vext- Treating a quantum mechanical system in a finite 
basis needs as a prerequisite a reasonable way to reduce the infinite dimensional Hilbert space 
to some subspace of finite dimension while keeping enough degrees of freedom to describe the 
essential physics. As shown above, an electron subjected to a magnetic field exhibits an 
equidistant energy spectrum of Landau-levels. The properties of a fractional quantum Hall 
ground state were shown to be reproducible when restricting the state to live in the lowest 
Landau level, neglecting admixture of higher Landau levels. Analytically this was done by 
Laughlin in From perturbation theory this limit corresponds to the case where the Landau 
level spacing ficoc is much bigger than the perturbation of the system, namely the Coulomb 
interaction and the impurity potential. 

This restriction to the lowest Landau level will also be applied here. But although it 
facilitates numerical calculations, some other problems are encountered when projecting to 
the lowest Landau level as described below in section 

The projection operator Plll can be constructed easily once the single-particle basis states 
|0, A;) = a|j|0) of the lowest Landau level from section [3.1.31 are known. Here denotes the 
creation operator of the single particle state ()19() with the momentum quantum number 
k G {0, . . . ,Ns — 1}. We can build many-particle states out of them as slater determinants. 
Alternatively we can write the state in second quantization 

\n,...,jNj=al...a]jO). (50) 

The dimension of the basis is (j^") and we can number the states by ascending multi- indices 
ji < . . . < Jn^, where Ng is the number of electrons, Ng the number of single particle states 
in the lowest Landau level. The projection operator can be constructed simply as 

Plll = \h, ■ ■ ■ ,jNe){ji, ■ ■ ■ JnJ- (51) 

jl<...<jiVe 
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The Hamiltonian from Equ. can now be projected by this operator. This yields 

n ( \ 

PlllHPlll = -^N^Plll + Plll + Yl " ^^ll- (52) 

\i=i i<j J 

Due to the projectors left an right to the operators, this Hamiltonian only causes mixing in 
the subspace of the lowest Landau level. Therefore, this is equivalent to diagonalizing the 
unprojected Hamiltonian in the many particle basis (|Sn|) . Stated differently, we just have to 
write the operators Vext and Vint in second quantization where the single- and two-particle 
matrix elements are calculated in the basis |0,j). The kinetic energy of ^"^^"^ is just a 
constant and does not affect the eigenstates at all. Omitting it is equivalent to changing the 
origin of the energy scale. This also makes clear why we have to diagonalize at all, since all 
many-particle states are degenerate with respect to the kinetic part of the Hamiltonian and 
the interaction cannot be treated as a small perturbation. To arrive at the eigenstates of this 
projected Hamiltonian we therefore have to diagonalize 



Hlll = ^ ^extji,j2'^iaj + ^ 'j^intji,j2,j3jiO-j^o-j2^jiO-h- (53) 

Jl>j2 = l il,j2,j3j'4 

Actually, there is an additional diagonal single particle operator in the case of Coulomb 
interaction. It is the constant interaction energy of one electron with its images in neighboring 
cells |H3 . The matrix elements for the interaction and the various external potentials used in 
the inhomogeneous system are found in the appendix IbI 

However, in this section we will first work with the homogeneous system, thus V^xt = 0. 

4.1.5 Groundstate energies and chemical potential 

The purpose of these calculations is to model a constriction in a fractional quantum Hall 
system and investigate how quasiholes near this barrier behave. More concrete, the question 
whether single quasiholes can pass through the barrier should be answered. The dependence 
on the barrier's parameters will be surveyed. To get a measure of the potentials needed to 
create a tunneling barrier (for electrons), the chemical potential of the electrons in the system 
must be calculated. If the height of the constriction is higher than the chemical potential, 
the electrons can only pass it by tunneling. The chemical potential is calculated for different 
system sizes from 4 to 6 electrons at a filling factor ^. Another purpose of this calculation 
is to confirm the system's incompressibility, which is recognizable as a jump in the chemical 
potential when going from slightly lower to slightly higher filling factors than ^. Tabled 
summarizes the ground state energies, their degeneracies and the gap of the system for both, 
Coulomb and hard-core interaction for different filling factors around i/ = |. Note that all 
energies are given in units of 

enu := , , , (54) 

where Ns is the respective number of flux quanta of the system. This energy unit will be used 
throughout the document. The energy per electron of the systems with N^/Ng =4/12, 5/15 
and 6/18 with Coulomb interaction are in agreement with the values calculated by Yoshioka 
in EH. 
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050057 


11 


205743 


4/13 H 


0.0000 


0.0 


13 


0.192552 


5/12 H 


0.609253 


0.121851 


12 


0.176295 


3/12 H 


0.0000 


0.0 


> 30 


n.c. 


5/15 H 


0.0000 


0.0 


3 


0.20662 


5/14 H 


0.21559 


0.043118 


14 


0.198772 


5/16 H 


0.0000 


0.0 


16 


0.199536 


6/15 H 


0.449654 


0.074942 


5 


0.235736 


4/15 H 


0.0000 


0.0 


> 30 


n.c. 


6/18 H 


0.000 


0.0 


3 


0.195 


6/17 H 


0.215296 


0.035883 


17 


0.133884 


6/19 H 


0.000 


0.0 


19 


0.197 


7/18 H 


0.5257 


0.0751 


18 


0.1304 


5/18 H 


0.0000 


0.0 


> 30 


n.c. 



Table 1: Ground state energies, degeneracies and gap of a system of 4 ... 6 electrons for a 
filling factor near v = The first column specifies the number of electrons, the 
number of flux quanta and the kind of interaction potential (H for hard-core, C for 
Coulomb). 
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From the values in tabled the chemical potential of electrons could be calculated just by 
taking the difference between ground state energies of a system with and electron. 
However, this approach has some drawbacks which are partly due to the small size of the 
system, partly because of the other quantities that are affected by changing N^. If we take 
a system with only four or five electrons and change the number of electrons by ±1 this 
change can hardly be treated as infinitesimal with respect to the total number of electrons, as 
demanded by a reasonable definition of /x. In fact, in the case of ^ — > ^ changing the number 
of electrons by +1 results in a filling of |, which should itself show a fractional quantum Hall 
effect. The second reason not to calculate the chemical potential this way is of more practical 
nature. When calculating the Coulomb interaction between electrons, to conserve the overall 
charge neutrality, the interaction with a uniform positive background of density cj+ = 
is taken into account. Thus, changing A'e also alters this background charge and affects the 
eigenenergy. 

To overcome the problems mentioned above, Yoshioka uses in [HSl a different approach to 
calculate the chemical potential for electrons in his numerical finite size studies by tracing it 
back to the dependence of the ground state energy on the number of flux quanta. To show 
the step in the chemical potential, around the FQH-effect fillings, the chemical potential is 
calculated for slightly lower (— ) and slightly higher (+) fillings of the system. This /i^ is 
calculated in the following way. Let Eq{u) be the energy per electron at the filling factor 

f^electron " ^GS ~ ^GS \^^) 

= {Ne + l)EQ{v{Ne + l))-N,EQ{u) 

BE 1 

= {N, + l)EQ{u) + N,^—-NeE^{v) 

JMs Of 

Eo{jy+)-Eo 

~ Eo[i^) + 1^ — ^ — 

h'^ — h> 

analogously (56) 

Eo{iy-)-Eo 



^'electron ~ + ^ 



— V 



Ne iV. N, 

where u = — — 1/4 



Ns ^ iV. -1 iV. + l 

Since here only the number of flux quanta is varied, the problem of changing the neutralizing 
background is not present. This calculation yields the values from table |2l for the chemical 
potential for adding an electron at slightly higher and slightly lower fillings than |. To address 
the calculation for hard core interaction, first note that according to section l4.1.31 the ground 
state energy in the system with hard-core interaction vanishes whenever the wavefunction has 
a threefold zero if two electrons come close to each other. Therefore, at least mNg zeros or 
flux quanta are needed. If there are less, we arrive at a finite ground state energy. On the 
other hand, if there is one flux quantum more than needed to establish these correlations, the 
hard-core interaction vanishes, too. So, adding a flux quantum is a gapless excitation of the 
system with hard-core interaction, while removing one needs a finite creation energy, since 
the correlations have to be changed. 

In the case of Coulomb interaction, there is a finite but yet small amount of energy to pay 
in either case. The negative chemical potential fi- means that the system wants to absorb 
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System 


f^electron ido{Ns=l)^ 


l^electron (elo{Ns=l)^ 


4/12 C 


0.27732 


-1.1745891 


5/15 C 


1.2132899 


-2.091855 


6/18 C 


2.1070906 


-2.4859929 


4/12 H 


7.304929 


0.0 


5/15 H 


11.277894 


0.0 


6/18 H 


15.528201 


0.0 



Table 2: Chemical potential for electrons and quasiholes for Coulomb (C) and hard core (H) 
interaction 

electrons until it has reached a filling of |, when the chemical potential jumps to a positive 
value and stops this process. 

Since the elementary excitations of the system are quasiparticles according to Laughlin 
we can relate this jump in the chemical potential to the creation energy of quasiparticles. 
Increasing the density amounts to adding quasielectrons, decreasing it can be understood as 
adding quasiholes. Since the quasiparticles only have ^ of an electron's charge, the creation 
of a quasihole needs eh = — ^/U_, for a quasielectron we need eg = Therefore, for both 

interactions we can conclude that the system is incompressible, because enlarging the area is 
equivalent to inserting quasiholes, while compressing it amounts to inserting quasielectrons, 
which in both cases costs a finite amount of energy for infinitesimal changes of the density. 

4.1.6 Correlations in homogeneous systems 

For electrons interacting with each other by the hard-core pseudopotential, the ground state 
should exhibit Laughlin- like correlations. This can not only be confirmed by vanishing hard- 
core interaction but also by looking at the two particle correlation function 

9(z) = lv^^(^lE'^(^"'+^^-^^-)l^)- (57) 

This correlation or pair distribution function was already used by Yoshioka in Since 
we are calculating in a basis with quasi-periodic boundary conditions (see section r3.1.4() this 
operator's matrix elements are a bit different from those used by Yoshioka. They are given 
in the appendix IbI 

The Laughlin wavefunction has a relative factor of (zj — Zj)"^ for every pair of electrons. In 
the case of = i = — it should lead to a correlation function q(z) that behaves like IzP'". 
Evaluating g{z) for each of the 3 degenerate ground states in a system with 5 electrons and 
15 flux quanta yields figure |SJ The right plot shows the quotient j^- The finite value it 
approaches for small z directly confirms the Laughlin- like correlations. 

4.2 Quasihole excitations 

A quasihole excitation is an excitation of the fractional quantum Hall system to which we 
can attribute particle like properties such as charge and position. The charge of this object 
was shown to be a fraction of the charge of an electron by Laughlin in Another argument 
than the plasma analogy used by Laughlin will be invoked here to make this fractional charge 
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g(z) 



mm" 




Figure 5: Two particle correlation function for a system with N^/Ns = 5/15 and hard-core 
interaction. The left plot shows the correlation function itself, the right plot is the 
quotient g{z)/\z\^. The finite value for \z\ — > confirms that g{z) oc \z\^ for small 
\z\ as expected for Laughlin-like correlations. 
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Figure 6: Gedankenexperiment to create a quasihole: A solenoid is pierced through the plane 
of the 2DEG in the gapped groundstate at filling = ^ and its magnetic flux 
is increased adiabatically by one flux quantum <I>o- The induced electric field E 
doesn't cause transitions but the charge of — ez/ is repulsed out of the area A. 



more plausible and to get a coarse idea and intuitive understanding of the connection of flux 
quanta and quasiholes. It is the following picture adapted from Shankar j30j . 

Imagine electrons confined in an infinite plane with a perpendicular homogeneous magnetic 
field B. Let the filling factor be the fraction u = —, m being an odd integer. The system 

2 

is then in a fractional quantum Hall state and shows a Hall conductivity of = while 
the longitudinal conductivity is zero, axx = 0. The ground state of the system is separated 
by a gap from the excited states. This gap prevents the system from undergoing transitions 
in the presence of an adiabatic perturbation. A thin solenoid is pierced through the surface 
as depicted in Fig. IHl 



The flux $ through this solenoid — initially being — will be increased adiabatically to one 
flux quantum $0 = 7- Since this is an adiabatic change, the system's fractional state will not 
be destroyed because it is gapped and the conductivity remains as stated above. The time- 
dependent change of ^ produces an induced field E by Maxwell's equation V x E = —dtB. 
This field is directed tangential to a circular boundary dA enclosing the solenoid (see fig. 
Due to the form of the conductivity tensor in the fractional quantum Hall state, this electric 
field produces a current flowing radially away from the solenoid through the boundary dA. 
This current can be integrated over time to yield the amount of charge Qrep that is repulsed 
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by the insertion of the flux quantum. This calculation yields 



rO 

Qrep = I dt I ds j (58) 



oo 



OA 




The total charge that flows through the boundary while increasing the flux in the solenoid 
in this process is Qrep — so there remains a positive charge deficit of Qqh — i^g inside the 
boundary. This is the charge that is attributed to the quasihole. As Shankar points out PU] . 
thinking in this picture, the fractional charge of a quasihole is a consequence of the quantized 
conductivity of the system and not vice versa. 



4.2.1 Variational wavefunction for a quasihole excitation 

Following the paper of Haldane and Rezayi we can construct a wavefunction that exhibits 
a quasihole excitation. Starting from a fractional quantum Hall system's ground state, a 
quasihole can be regarded as a zero in the wavefunction relative to all electrons' coordinates. 
As already stated in section l4.1.21 to fix a zero of the trial wavefunction at a specified position 
with respect to every electron, an additional flux quantum must be introduced into the system 
to gain the desired freedom and in accordance with the Aharanov-Bohm effect. 

Such a wavefunction of a state with a localized quasihole at zq = xq + iy^ was given along 
with the system's ground state in this paper [^. It can be understood as a modification of 
the ground state wavefunction (|45() . 

/ —i{z* — z*) L \ 

j=l 2 2 

j<k ^ ^ 

X F^.^i(Z + -)exp(^ j 

If we compare the ground state wavefunction (|45() with (|59|) . the procedure to create the 
quasihole excitation can formally be seen as a multiplication with a relative zero (since ??(z) ~ 
z for small \z\) for every electrons' coordinate Zi with respect to the location of the hole zq. 
Additionally, the center of mass function Fn^-^i{Z) satisfies the boundary conditions for A'^g + l 
flux quanta and is translated by ^ . In the Gaussian the number of flux quanta Ng is replaced 
by A'^s + 1. Physically this can be interpreted as inserting an additional flux quantum at zq, 
manifest in an Aharanov-Bohm phase of 2tt when encircling zq with one of the electrons. The 
increase of the number of flux quanta is also reflected in the change of the number of flux 
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quanta in the Gaussian. Finally the center of mass coordinate has to be translated which will 
appear to be crucial for conserving the boundary conditions. 

These operations will be incorporated into a quasihole creation operator Onoieizo). It 
will be derived such that acting on the ground state function results in the quasihole 
wavefunction (j59|) . 



After deriving this operator we can let it act on a many particle basis state of a system with Ng 
flux quanta. These many particle basis states are linear combinations of slater determinants 
of the single particle basis states from Equ. ()2U() in section 13.1.41 Since we are adding a flux 
quantum by this process, the resulting vector can be expressed in the basis of a system with 
A'^s + 1 flux quanta. If we once know how these basis states of the Ns system translate to 
those of the A^^ + 1 system with an inserted quasihole, the operator can as well be applied 
to arbitrary states expressed in the 'W<j"-basis. From the operator's definition it is clear 
that it always produces only trial wavefunctions for quasihole states which may be good or 
not. For homogeneous systems with short-range interaction or with Coulomb interaction this 
approach will be shown to produce excitations that are identical or close to the true quasihole 
excitation, respectively. 

The rather lengthy calculation of the operator's derivation can be found in Appendix 1X1 

4.2.2 Creating a quasihole by pinning of a vortex 

Regarding the trial wavefunction for a quasihole excitation in Equ. (|59)) from the previous 
section, it is possible to make two points. First, due to the m— fold zero in the relative 
coordinates this wavefunction must have a vanishing expectation value of the short-range 
interaction energy as was shown in section 14.1.31 Second, because of this vanishing ground 
state energy, all states only differing in the quasihole's position zq will be degenerate. Since 
hard-core interaction is a positive definite pseudo-potential (in spin polarized systems), the 
ground state of a system with Ng = mNe + 1 flux quanta cannot have lower energy than zero. 
Thus, Equ. must be a vector that lies in the space of the degenerate ground state of a 
homogeneous system. So, it must be possible to find a ground state with vanishing energy, 
even if the position of one vortex is fixed as a quasihole at zq. 

This constraint can be imposed to the system by adding a delta potential J2^=i ^i^i ~ 
xo)^{yi ~ yo) at Zq = XQ + iyo to the Hamiltonian. In the case of a filling factor of = |, diag- 
onalizing the Hamiltonian with short-range interaction H41() must then result in wavefunction 
(|59|) . since this is the only wavefunction that at the same time suffices all constraints: The 
periodic boundary conditions, a zero at zq and 3 relative zeros on every electron. The part 
-F/Vs+i of the wavefunction l59l still is defined by a wavevector K and 3 zeros Zi, (see eq. 
and hence produces a threefold degeneracy as in the homogeneous case. 

4.2.3 Quasiholes in homogeneous systems with hard-core interaction 

In this section quasihole excitations will be investigated for a homogeneous system of electrons 
interacting via hard-core interaction. The filling factor is set to u = ^, realized by Ag = 5 
electrons and Ns = 15 flux quanta. 

Applying the creation operator for a quasihole, deduced in section to each of the three 
degenerate ground states of the N^/Ng = 5/15 system results in three quasihole states having 




(60) 
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the hole at the same position. The electronic density of one of these states is shown in Fig. 
[71 left plot. The states are not identical and only after superposition of the three densities a 
circular symmetric result is achieved (apart from finite size effects, see below), as shown in 
Fig. [3 right plot. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 



Figure 7: Electronic density of a quasihole-state created by application of the quasihole cre- 
ation operator on the ground state of a homogeneous system with N^/Ns = 5/15 
and hard-core interaction. Left plot: Density after application on one of the three 
deg. ground states; Right plot: Densitiy obtained by superposition of the three 
degenerate states. 



To get a measure of the charge qgh repulsed by inserting the quasihole, a Gaussian of width 
Iq is fitted to the density. The area below this Gaussian is 2ttIq which takes exactly of 
the system's area. In a homogeneous system, this Gaussian would therefore envelop | of 
an electron's charge. Fig. IS] shows a quite good agreement of the the Gaussian (5 + ^)(1 — 

exp(— J^)) with the density profile which is an indicator of a fractional amount of charge 

</ : = ie that was repulsed by insertion of the vortex (compare section 14.2(1 . The slight 
asymmetry encountered in the density profile is probably a finite size effect. In an infinite 
system the cut must be symmetric. Still unclear is the stability of this excitation. The 
projection of the quasihole-state to the eigenstates of the homogeneous N^/Ng = 5/16-system 
gives the answer: The excitation lies completely inside the space of the 16-fold degenerate 
ground state (see figure hence it is an eigenstate itself and a gapless excitation as in the 
case of the addition of a flux quantum in section 14.1.51 By application of the quasihole creation 
operator on each of the three degenerate ground states of the homogeneous N^/Ng = 5/15- 
system, three linearly independent degenerate states with eigenenergy zero and a quasihole 
at the desired position could be created. Stated differently, the space of quasihole excitations 
is a three-dimensional subspace of the ground state sector of the N(,/Ns = 5/16-system. This 
is true for every point in the unit cell the quasihole is created. The vanishing eigenenergies 
in systems with hard-core interaction are an indicator for Laughlin-like correlations between 
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Figure 8: Cut along x = y through the electronic density of the N^/Ns = 5/16 system with 
inserted quasihole. A Gaussian of width Iq covers the area of one flux quantum and 
fits quite well the profile. For asymmetry see text. 
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Figure 9: Absolute values of projection coefficients of a quasihole state created by application 
of the quasihole creation operator on one of the ground states of the N^/Ng = 5/15- 
system by projecting on the eigenstates of a homogeneous N^/Ns = 5/16-system. 
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the electrons. Bearing the construction of the quasihole operator in mind, it is clear that 
these correlations — established in the ground state of the homogeneous = | system — are 
conveyed by the operator to the state containing a quasihole. 

As described in section 14.2.21 there is an alternative way to insert a quasiholes by pinning 
a flux quantum through a delta-shaped potential at a fixed position. In the absence of 
other external potentials and using hard-core interaction this results in a threefold degenerate 
ground state with vanishing energy, as expected. The density of one of these states looks very 
similar to the result gained earlier by means of the operator (see left plot in Fig. CJ. Projecting 
this quasihole state to the eigenstates of the homogeneous N^/Ng = 5/16-system shows that 
this excitation lies completely inside the space spanned by the 16-fold degenerate ground state 
vectors (figure ITU)) . Therefore the stability of this excitation is obvious since it is an eigenstate 
again. This is of course independent of the position where the hole was created. 
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Figure 10: Absolute values of projection coefficients of the quasihole-state created by diago- 
nalizing a N^/Ns = 5/16-system with hard-core interaction and a delta potential 
by projecting on the eigenstates of a homogeneous N^/Ns = 5/16-system. 



Another interesting feature of the method using a delta potential to create quasihole excita- 
tions are the excited states of such a system. As already known from Tab. Q in a 5/16-system 
there is a 16-fold degenerate ground state. Placing a delta potential into such a system, which 
is chosen to have a much smaller integral energy than the gap of 0.200enu, should cause a 
mixing mainly within the ground state sector. Compared to the gap of a homogeneous 5/15- 
system, the gap of this system is almost unaltered, which confirms that there is negligible 
mixing between states below and above the gap. The diagonalization with a delta-potential 
^delta ~ 0.01enuLiL2 X^^''^ 6{ri — fo) at the origin (fo = 0) yields the spectrum in Fig. ^2 
The threefold degenerate ground state with zero energy is confirmed. The 16-fold ground 
state splits into three different energy bands, two of which are made up out of three, one out 
of ten almost degenerate states. Evaluating the densities for these states reveals their nature: 
The 10-fold quasi-degenerate subspace consists partly of states looking like superpositions of 
quasiholes located at positions that are different from tq. The best example is state 11, whose 
density is given in Fig. El (right plot). The quasihole seems to be localized mainly at the 
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origin. The shght dip in the density at rg is a common feature of these 10 states, but most 
of the other 9 eigenstates have a more comphcated density profile. Nevertheless, another 
sign that supports the interpretation of these states to be quasihole ground states at different 
positions will be given in section 

A qualitatively different picture is obtained for the three higher states 13,14,15. The super- 
position of their densities shows a ring like structure around the delta's position, as found in 
Fig. ^l(left). This looks very much like an excited state of the quasihole. The inner radius 
of this ring is about 21q, the outer one about 3/o- 



spectrum 5/16 with delta 
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Figure 11: Lowest 20 eigenenergies of the 5/16-system with hardcore interaction and a delta 
potential at rb = 0. The gap of a homogeneous 5/15-system is given for compari- 
son. 
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Figure 12: Left: Superposition of the densities of the threefold quasi-degenerate states 
13,14,15 in a 5/16-system with a delta potential at vq = (0.5,0.5). Right: Density 
of eigenstate 11 of the same system. The magnetic length is about /q — -jq- 
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4.2.4 Quasiholes in homogeneous systems with Coulomb interaction 



Again a system of 5 electrons is chosen, the magnetic field causes a flux of 15 flux quanta 
through the unit cell, but the electron interact via the Coulomb potential. The same quantities 
as in the case of hard-core interaction are investigated. Beginning with the density, an obvious 
difference between the quasihole created by delta potentials in contrast to those obtained by 
annlication of the onerator ('comnare Fiff. 11 .'-{1 ('left') arid Fiff. 1141 1. is. that the former ones 

' I 1 T 1 1 1 I I ' I 1 1 1 ] 1 I — I 





Figure 13: Density of a qu'dsihole in a 5/16-system with Coulomb interacfen created by appli- 
cation of the creation operator to a ground state of the homogeneous 5/15-system. 
Left plot: One single state; Right: Superposition of the densities of three states 
obtained from the tree deg. ground states. 
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Figure 14: Density of a quasihole in a 5/16-system with Coulomb interaction and a delta 
potential obtained as the ground state by diagonalization. 



quasihole is by looking at their spectral decomposition. The states are projected onto the 
eigenstates of a homogeneous 5/16-system. This reveals another difference: In the case of 
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Figure 15: Absolute values of projection coefficients for a quasihole-state with Coulomb inter- 
action obtained by projecting on the eigenstates of a homogeneous 5/16-system. 
Left plot: Quasihole is generated by the creation operator. The excitation energy 
is 0.0023975enu above the GS. Right plot: Quasihole obtained as ground state 
from diagonalization with a delta potential. Its excitation energy is 0.00079831enu 
above the GS. 



Coulomb interaction the quasihole generated by a delta potential has smaller contributions of 
the excited states than the one generated by the operator. This can be attributed to the fact 
that the operator was derived from Laughlin's trial wavefunction, which is exact for hard-core 
interaction, but which in turn is only a good approximation in the case of Coulomb interaction. 
Aside from the projection coefficients \Ci\ themselves, the energy of the quasihole states can 
be computed by using Equ. H61|) and the known eigenenergies Ei of the homogeneous system 

E = Y,\Ci\^Ei. (61) 

i 

Tab. Olshows the results for both — Coulomb and hard-core interaction. Obviously, quasihole 
states are no longer eigenstates of the system when Coulomb interaction is in use. Further- 
more, the delta potential created quasiholes have excitation energies which are lower by 
approximately a factor of three compared to those obtained by means of the operator. The 



System 


Method 




5/15 C 


operator 


0.0023975 


5/15 C 


delta potential 


0.0007983 


5/15 H 


operator 


0.0 


5/15 H 


delta potential 


0.0 



Table 3: Comparison of the excitation energy (above GS) of a quasihole created by application 
of the operator to a homogeneous 5/15-systems's ground state and by diagonalizing a 
5/16-system with a delta potential for hard-core (H) and for Coulomb (C) interaction. 

stability of a hole excitation in the system with Coulomb interaction is not obvious in ad- 
vance. To clarify this, the time evolution of the occupation numbers nj{t) is surveyed. This 
quantity can be calculated according to Equ. (|62|) . once the spectral decomposition of the 
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quasihole-state into eigenstates of the homogeneous 5/16-system and thus its time evolution 
^(t) is known, 

n,{t) = {^it)\alaj\^{t)). (62) 
Fig. ^1 shows the time dependence of nj{t). The color code indicates the occupation 




quantum number k 

Figure 16: Time evolution of the occupation numbers nj{t) for the quasihole obtained by 
diagonalization with a delta potential at tq = (0-5, 0.5) and Coulomb interaction. 

probability of the single particle state j (x-axis) at time t (y-axis). In a homogeneous 5/15- 
system the occupation of a single particle state is equal |. In contrast in Fig. ^Jthe states 
around j = 6 . . . 9 have a lower occupation probability due to the inserted quasihole. This dip 
in the occupation is nearly constant in time and only shows some minor oscillations. Since 
the quantum number j is coupled to the point Xj = around which the single electron 

state is localized, we can conclude from the Constance of nj{t) that the hole does not move in 
x-direction. As both axes - x and y - of the unit cell in this particular system differ only by the 
selected gauge, we can infer that there can neither be a motion of the quasihole in y-direction. 
Instead it rests in the middle of the cell where it was created and shows some breathing mode 
like oscillations. These can be observed by calculating the density as a function of time. 

4.3 Conclusions on homogeneous systems 

Laughlin's trial wavefunctions for the ground state and quasihole-excited state in a homo- 
geneous system at = | were introduced and shown to be eigenfunctions for a special 
short-range interaction. 

In a homogeneous system known results were reproduced thereby verifying our calculations. 
This yields eigenenergies from which the chemical potential /i for electrons is calculated and 
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the incompressibihty of the system is confirmed for both — Coulomb and short-ranged — 
interaction. The correlations were found to be identical to those of Laughlin's wavefunction 
and the short-range interaction energy {Vgliort-Tange) GS is understood to be a suitable in- 
dicator for this. Based on the trial wavefunction of a quasihole excitation adapted to our 
geometry, a quasihole creation operator could be derived. Its validity was verified by look- 
ing at the properties of the excitation it produces upon application to the ground state of a 
homogeneous = | system: The excitation carries a fractional charge of and lies in the 
ground state sector (as expected) for short-ranged interaction, thus having excitation energy 
0. For Coulomb interaction a finite excitation energy is found. 

An alternative method of creating a quasihole is by fixing one vortex of the many-particle 
wavefunction in a system with one flux quantum in excess (iV<j = 3Nf. + 1); this was possible 
by diagonalizing this system with a delta potential to pin one zero of the wavefunction. 
Comparing the ground state obtained this way to the excitation created by the quasihole- 
operator reveals both approaches to be equivalent for short-ranged interaction. In the case 
of Coulomb interaction an energetically lower (factor of 4) quasihole excitation is obtained 
by the second procedure. Small contributions above the ground state of the homogeneous 
system to this state cause a breathing mode oscillation, but the quasihole is still stable. 

Another interesting feature of a delta potential in a system with Ng = SNg + 1 flux quanta 
are localized quasihole states at the potential's position. These will be of use in the following 
section to create a setup where tunneling can be observed. Furthermore, excited states of 
the quasihole at the delta potential with a ring-like structure were found. To some extent, 
results from section |21 indicate these states to be qualitatively different from the lower lying 
ones. Analyzing them more thoroughly could be the matter of further work. 
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5 Tunneling of quasiholes between delta potentials 

One aim of this work is to investigate the ability of quasiholes to tunnel through a constriction 
in the system. Before turning to this more complex problem, a much simpler case will be 
investigated here. The question is whether there is evidence for quasihole tunneling in our 
system at all. A setup that is helpful to shed light on this was already found in section 
14.2.31 A system of A'^e electrons interacting via hard-core interaction, Ng = 3Ne + 1 flux 
quanta and a delta potential. The "additional" flux quantum was found to have a threefold 
degenerate localized state at the delta potential's position. Introducing into this system 
a second delta potential at a different position should result in tunneling of the quasihole 
between the localized states at the one delta and the other. Analogous to the picture of a 
single particle tunneling through a potential barrier, the role of the barrier is played here 
by the space between the delta potentials where the quasihole has to pay much energy to 
live. If tunneling occurs, symmetric or antisymmetric linear combinations made up out of 
the localized states should be observed to form the ground state and lowest excited states, 
respectively. 

The idea behind this construction is to assume the system of Ne electrons and Ng = SN^ + 1 
flux quanta to be effectively treatable as a single-quasihole system. One has however to 
remember that this single quasihole lives on the background of the homogeneous system's 
fractional quantum Hall state. This state can be thought of as the vacuum with respect to 
quasiholes. Returning to the single quasiparticle picture, let us denote the l-th localized state 
of the quasihole at the delta potential i as \l)i, where i = 1,2 and I = 1,2, 3. Of course, in a 
system with two delta potentials these states are no longer eigenstates. Instead we use them as 
a basis and try to obtain a reasonable estimate for the ground state and lowest excited states 
within this subspace. The restriction to this space can be justified by perturbation theory, 
since due to the symmetry of the problem with respect to exchanging the delta potentials, 
the states \l)i must be degenerate. They have some energy e = {l\iH\l)i for i = 1,2 and 
I = 1,2,3, H being the Hamiltonian of the system with both delta peaks. The mixing due to 
tunneling must therefore be most dominant for these states and we will take it into account 
by a Bardeen-type tunneling Hamiltonian Ht = ;/ i/,/'|/)i2(^'| + ^^/|^')2l(^|• 
Now, there needs to be an assumption to set up the tunneling coefficients ti^i'. To render 
a connection between the quantum number / and some physical observable, one has to keep 
in mind that the threefold degeneracy of the quasihole state emerged from the (threefold) 
degeneracy of the ground state of the homogeneous system aX u = ^. The latter degeneracy 
originated from conservation of the center-of-mass momentum Py. Returning to the system 
with a quasihole and one delta potential, there are three orthogonal quasihole states at every 
point zq in the unit cell. So there must be another operator — call it G{zq) (like the cen- 
ter of mass Py in the homogeneous case) — that commutes with H and whose eigenstates 
discriminate three orthogonal eigenvectors within this subspace. If G is independent of the 
point Zq, we will be able to separate the quasihole state sector into three eigenspaces of G 
and label them by I = 1,2,3. Of course, it is possible that G{zq) depends on the position 
Zq of the quasihole. But however it is tempting to assume that it does not. Then it would 
be reasonable to suppose that the tunneling coefficients t/^// only depend on whether the two 
coupled states \l)i and \l')2 lie in the same eigenspace of G, / = or if they lie in different 
ones / 7^ This assumption yields tunneling coefficients ti^i/ = 5i^i't -|- (1 — 6i^i')ti t,ti G C 
The Hamiltonian H = e + Ht can be written in the basis b of the direct sum of the ground 
state spaces of a quasihole at the one delta or the other: b = {|l)i, |2)i, |3)i, |1)2, |2)2, |3)2}, 
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where the states are labeled by the eigenstates of the suppositional operator G. Hence the 
Hamiltonian has the following matrix form 



Its eigenspectrum is 
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If |t + ti| > then A > K and the spectrum has the structure consisting of one lowest 
eigenstate e — A followed by a two-fold eigenstate e — another two-fold e + k and a highest 
one with e -|- A. This structure will be recognized in the numerically obtained spectrum in the 
next section. 

5.1 System with two delta potentials and one flux quantum in excess 

In a system of 5 electrons and 16 flux quanta there is a 16-fold degenerate ground state in 
the homogeneous system, as seen from Tab. ^ Two delta potentials shall be introduced into 
this system. Their integral potential energy is chosen very small, such that they will not 
cause strong mixing between states below and above the gap. Introducing only one delta of 
height 0.01LiL2enu produces the spectrum shown in Fig. ^Jof section l4.2.Hl The splitting of 
the 16-fold degenerate ground state is below half of the gap's size. If we are to introduce the 
second delta and want the mixing above the gap to be negligible, the amplitude of both deltas 
has to be yet a bit smaller to be on the save side. Here we chose a height of 0.005ennLiL2 
for both deltas at the positions (0.0, 0.5) and (0.5, 0.5) in the unit cell. Diagonalization yields 
the spectrum seen in Fig. El red plot. The ground state sector exhibits a splitting, where 
the lowest six energies form a band that is supposed to originate from the coupling of the 
two times three degenerate ground states for the quasiholes localized at the two deltas, as 
described before. These states should therefore be symmetric or antisymmetric superpositions 
of eigenstates for the quasihole at the left or the right delta. The structure of this spectrum 
is the same as expected from the simple model (Equ. I65() . A single ground state, two pairs 
of two fold degenerate states and a single highest state. The density of the ground state is 
given in Fig. 1181 It shows two dips at the positions of the delta potentials which in turn 
are located in a "trench" of low density connecting the two potential peaks. This picture 
suggests, that at this distance (~ 5/o) between the delta potentials, the delta peaks still can 
not be treated as decoupled. Therefore a splitting of the degenerate states is already visible at 
this distance, as seen from the spectrum in Fig. 1171 red plot. To determine if this state \^) is 
a superposition of a quasihole ground state located left and one located right, the projection 
coefficients of l^*) onto the eigenstates of a system with one delta at (0.0, 0.5) (called system 1 
in the following) and onto those of a system with a delta potential at (0.5, 0.5) (called system 
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Figure 17: Spectrum of a 5/16-system with two delta potentials 0.5Li apart; red plot: Both 
delta potentials have a prefactor of 0.005LiL2enit; green plot: The prefactor is 
0.004 and 0.00QLiL2enu respectively. 
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Figure 18: Density of a 5/16-system with two delta potentials (strength 0.005LiL2enu) being 
0.5Li apart. The two dips in the density (location of the potentials) are connected 
by a "trench" of low density. 
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Figure 19: Projection coefficients for the state with two delta potentials onto the eigenstates 
of the system with only one delta potential at the left or the right position. The 
distance between the potentials is 0.5Li. The gap of the homogeneous system is 
0.2enu. 



2 in the following) are calculated. Due to real vectors (for this special choice of positions), the 
coefficients in Fig. EJare real, too. A first point to make is that contributions of states above 
the gap (0.2enu) are small. If the ground state |^) of the system with both deltas would be 
a symmetric or antisymmetric linear combination of the ground states of the systems with 
a single delta and if the eigenspaces of localized quasihole states at different positions were 
orthogonal, one would expect half of the state to lie in the ground state sector of system 1, 
the other half to lie in the ground state space of system 2. Unfortunately these spaces are not 
orthogonal. This non-orthogonality results from the fact that a quasihole excitation created at 
an arbitrary positions in the unit cell always lies in the 16-fold degenerate ground state space 
of the homogeneous system. Due to the symmetry upon exchanging the two delta potentials, 
the projection coefficients however must have the same magnitude for system 1 and system 2 
(if the corresponding state on which is projected is non-degenerate). To take into account the 
overlap between the threefold degenerate ground state spaces of system 1 and system 2, the 
amount of the wavefunction that resides in the product space can be calculated as follows. 
First the projection of |^) onto the eigenspace of system 1 is performed. The result is given 
in Fig. The part of the wavefunction that is projected by Pgsi into the ground state 
sector will be called I^'gs'i) = -Pgs'iI^)- It makes up 60.5 per cent of the state. The "rest" of 
the state not lying within this space, (1 — Pgsi)^, is of course orthogonal to the ground state 
space 1. It will be projected in a second step onto the eigenspace of system 2. The resulting 
absolute values are given in Fig. 12111 A big part \"^gS2±gsi) = PgS2{^ — Pgsi)]"^) lies in the 
ground state space 2. It amounts to 34.5 per cent of the whole state. The total amount of 
1^) lying in the sector spanned by both ground state eigenspaces of system 1 and system 2, 
respectively, is just given by the squared norm iV^ = (^-gsiI^gsi) + {'^gS2±gsi\^GS2±gsi) 
and results in 94.4 per cent. Thus the assumption for the state to be mainly an equally 
weighted linear combination of the ground states of the quasihole localized at the left or at 
the right delta potential is confirmed. 

Another interesting finding is that there are nearly no contributions from the states 14, 15, 16 
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Figure 20: The contributions of excited states of system 1 to the ground state projected onto 
the eigenstates of system 2. A high overlap with the ground state is found. The 
states 14,15,16 around Ei — Eqs — 0.76 do not hardly contribute. 

around — Eqs = 0.76enu in Fig. 1191 This is another indication that these three states are 
qualitatively different from the ten states around Ei — Eqs — 0.035enu, as already proposed 
in section 14.2.31 

To determine whether the ground state is a symmetric or antisymmetric linear combina- 
tion of the ground states of system 1 and system 2, looking at the phases of the projection 
coefficients cannot give an answer to this question. The reason is simply that the phase of an 
eigenvector in system 1 is arbitrary. The phase — or the sign in our case of real vectors 
— of the coefficients in Fig. EJthus depends on the arbitrary choice of the respective eigen- 
vector's phase. The only meaningful definition of a symmetric or antisymmetric wavefunction 
must take into account the form of the contributing eigenstates of system 1 and system 2 in 
x-space. It has to be checked if the superposition is an odd or an even function with respect 
to interchanging the delta potentials. For this purpose Fig. 1181 which is the electronic density 
can be regarded as the density of the quasihole wavefunction. To get a quasihole-density from 
this plot, one has to subtract from the mean electronic density (5 in the units used here) this 
profile: The minima in Fig. 1181 become maxima and vice versa (inversion of the gray scale). 
What was called "trench" earlier now becomes a "ridge". This shape of the wavefunction 
without a node between its two maxima one would call a symmetric linear combination or 
a binding orbital in atomic physics. This qualitative shape of the quasihole wavefunction 
is more easily recognizable in Fig. I2H where a cut through the density along the x-axis is 
shown. 

5.2 Distance dependence of the tunneling 

For smaller distances between the delta peaks the coupling due to tunneling is expected to 
increase because the width of the "barrier" decreases. This dependence will be investigated 
here. In Fig. [^a cut through the density parallel to the x-axis at y = O.5L2 is plotted for 
different positions xi of the second delta peak. The first one is kept fixed at (0.0,0.5). For 
distances > 0.26 the plots reproduce the two dips at the respective positions of the deltas. 
For smaller separations, there is a qualitative change and the curve for xi = O.lLi only shows 
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Figure 21: Cut through the density along the line between two delta potentials in a 5/16- 
system for various positions xi of the second delta. The first delta is fixed at 
xq = 0. A qualitatively change occurs below xi ~ 0.2Li ~ 21q where the two dips 
merge into one. 



one dip located in the middle of the two potentials. This qualitative change occurs at a 
separation of 2/o — 0.2Li. It is understandable, since the spacial extent of a quasihole is of 
the order Iq (width of the Gaussian fitted to the density, see Fig. IB}. If the delta potentials 
come closer than twice this extention the states' overlap of localized quasiholes left and right 
becomes very large. The state that is preferred then is a quasihole situated in between both 
delta potentials. This can be infered from the density Fig. [^as well as from the spectrum 
depending on the potentials' separation Fig. 1221 For x = the spectrum evolves into the 
state of a localized quasihole at the origin. For small separations between the potentials (i.e. 
< 21q) the energy of all states rises, but the structure of the spectrum remains the same. The 
rise in energy can simply be attributed to the potential energy caused by the delta potentials 
sitting on the slope of the quasihole's state where the wavefunction has small modulus but is 
not zero. At a separation around > 21q the symmetric or antisymmetric state of the tunneling 
begins to be favored. 

The dependence of the eigenspectrum on the distance between the two delta potentials, 
Fig. 1221 depicts at the point x = O.SLi the situation from Fig. 1171 red plot. The six lowest 
eigenenergies are very close to each other. Decreasing the distance increases the splitting 
between them. The comparison to the model at the beginning of this section, which lead to 
Eq. 1651 shows that the spectrum evolves a bit more complex than expected. The highest state 
does not split off as much with respect to the higher pair of degenerate states as the lowest does 
to the lower pair (as it should be according to the model). Therefore, the model made above 
is not applicable. This might be caused by mixing with higher states but the lowest six, which 
was neglected before. These higher states correspond - as stated above - either to quasiholes 
being localized at different positions than those of the two deltas under consideration or to 
excited states of the quasihole at the deltas' positions. Another reason might be the influence 
of the delta potentials on the "vacuum" state, i.e. the fractional quantum Hall state without 
any quasihole, which is the background on which the single quasiparticle lives. It could cause 
deviations via influencing the tunneling matrix elements. If there are more than the proposed 
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Figure 22: Spectrum of a 5/16-system depending on the spacing between two delta potentials 
(strength = 0.005LiL2enu each). At x = the two delta potentials effectively act 
as a single one. Left: lowest 7 eigenstates. Right: Splitting between the 5. exc. 
and the ground state. A qualitative change for x < 0.21q can be observed. 



two different coefficients t and ti, the calculated spectrum might be reproduced. Although 
the simple model in the beginning of this section does not apply unmodified, it is evident that 
the splitting of the six lowest states is due to the tunnel-coupling. The dependence of the 
energy difference between the highest of these states and the ground state energy produces 
Fig. 1221 The splitting increases with reduced distance between the two deltas until a distance 
of 0.2Li ~ 21q at which, as said before, the quasihole begins sitting in between the two delta 
potentials. This distance sets the lower limit for the tunneling regime on the spacing between 
the potentials. 

5.3 Asymmetric delta potentials 

Even at a distance of 5/o between the delta potentials there is still a coupling of the states 
with a localized quasihole on the left delta and those where the quasihole is pinned by the 
right one. To confirm, that this tunnel coupling is distance dependent without resorting to 
the spectrum, the following modification is made. The symmetry of the system is broken 
by increasing the height of one delta, while decreasing the other. The higher one, located 
at (0.0,0.0) will have an amplitude of 0.006-LiL2enii, the lower one at (0.5,0.0) has a height 
of 0.004LiL2ereM. These values were chosen such that the difference in the peak height is 
a bit larger than the splitting of the lowest six energies observed in Fig. 1171 red plot. The 
asymmetry in the "binding energy" of 0.002enu between a quasihole localized at the left and 
one at the right delta potential was chosen such that at a spacing of x = 0.3, the tunnel- 
splitting of about 0.045enu (Fig. 122(1 should dominate the process. The spectrum for the 
asymmetric case is given by the green points in Fig. Eland is very similar to the symmetric 
case, but the splitting between the six lowest states should now correspond to states where the 
quasihole is either localized at the left, stronger delta, giving lower energy or at the right one, 
resulting in a higher eigenenergy. A cut through the density (averaged over degenerate states 
2,3 and 4,5) can be found in Fig. 1231 The asymmetry originating from the different peak 
heights is quite pronounced. The three lower lying states 1,2,3 have the quasihole stronger 
localized at the larger delta potential at x = 0, while the three higher ones ~ 4,5,6 ~ show a 
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Figure 23: Cut through the density along the x-axis for a 5/16-system with two delta poten- 
tials, one is fixed at (0.0,0.0). Its strength is 0.006-LiL2enu. The second delta has 
a strength of 0.004LiL2enn. Left plot: The second delta is situated at (0.5,0.0). 
Right plot: The second delta sits at (0.3,0.0). 



stronger dip in the density at the delta at x = 0.5Li. Decreasing the distance between the 
deltas must increase the tunneling and therefore favor the symmetric or antisymmetric state 
if the tunnel-coupling exceeds the asymmetry in the deltas' peak heights. For a distance of 
X = 0.3Li this results in density profiles like in Fig. 1231 The asymmetry between the deltas 
is still strongly visible in the lower three states, while the upper three states are already very 
symmetric. A qualitative difference between the lower and the higher three states now can 
be attributed less to the place where the quasihole is localized, but to the symmetry of the 
states: The lower three ones favor a low density between the deltas, while the upper three 
ones have a peak in the density in between. This situation resembles the case of a coupled two 
level system, where we have symmetric and antisymmetric linear combinations of the single 
systems' eigenfunctions. But still the asymmetry in the density between both delta potentials 
is visible. In the limit of much stronger tunneling, the state must become symmetric. In the 
ideal case, from the calculation with asymmetric delta potentials, it is possible to extract 
information about the distance dependence of the tunneling matrix element t. Starting at 
large separations, tunneling is negligible and the quasihole state localized at the higher delta is 
favorable. Decreasing the distance until the tunneling becomes dominant over the asymmetry 
we should — at sufficient strong tunneling — end up in a symmetric density profile (as seen 
in the symmetric case). The reason why this cannot be seen clearly in our system is due to its 
finite size: Even at a separation of half a unit cell, the tunneling is still quite strong. This was 
already anticipated from the density profile in Fig. ^1 ( "trench" in the density between the 
potentials), from the finite splitting of the ground state sector in Fig. |221at x = 0.5 and finally 
also confirmed by the density profile in Fig. I23L If the tunneling was negligible, the latter 
density plot would have to show the quasihole localized at the higher delta for the ground 
state, while it should reside at the lower delta in the excited state. Although this tendency 
is visible, it is already superposed by the effect of tunneling. The energy splitting between 
these states must be comparable to the difference in the peak height, which is however true 
(compare Fig. 1171 green plot). This lower limit on the tunnel coupling given by the finite size 
of the system demands a stronger asymmetry in the peak height. On the other hand, going 
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to small distances between the potentials, the tunneling regime was seen to be limited at a 
separation of ~ 21q (saturation in Fig. 1221 localized quasihole between the deltas in Fig. 121(1 . 

5.4 Conclusions on tunneling of a quasihole between delta potentials 

Two delta potentials in a system with an additional flux quantum are found to be a good 
model system to get a view on tunneling of quasilioles. The ground state of this system could 
be shown to be to a high degree (94.4 per cent) a linear combination with equal coefficients 
of the quasihole ground states localized at the potentials, thus giving an evidence that the 
description in terms of a single-quasiparticle picture is valid. Although the spectrum becomes 
complicated for strong tunneling, for weak tunneling its lowest energies resemble the structure 
obtained by the simple model just considering tunneling between the lowest localized quasihole 
states. 

As a side effect, the classification of the excited states below the gap (see section 14. 2. 3(1 
could be substantiated: States with a ring- like structure were found to be absent in the 
spectral decomposition of the tunnel setup's groundstate, whereas all states with lower energy 
contributed. 

The expectation was confirmed that tunneling should increase with decreasing distance 
between the potentials. It would be interesting to compare this distance dependence with the 
overlap of the quasihole wavefunctions. The saturation of the tunnel splitting for distances 
smaller than 21q however tells us, that the dependence of the tunneling on the overlap can 
only hold for sufficient large separations. 

The finite size of the system was found to give a lower limit on the strength of the tunneling, 
whereas the finite spacial extent of the quasihole wavefunction on the other side dictates an 
upper bound for it. These limitations prevented the calculations with asymmetric potentials 
to yield quantitatively usable results. However, these limitations were not exhausted yet. As 
seen from Fig. ^Jin section (4.2. 3! there is an excited eigenstate of a quasihole localized mainly 
at the position diametrically opposed to the delta potential. Owing to its orthogonality 
to the quasihole ground state, choosing two diametrically opposed points within the unit 
cell to place the delta potentials will probably result in weakest possible tunneling (as also 
intuitively expected). Therefore one could start with a smaller asymmetry in the peak height 
and possibly extract the distance dependence of t from the point of crossover between the 
asymmetry-dominated and the tunneling regime. 

Possible further work could also include a comparison of the results obtained here to those 
for electrons (or holes). The same calculations at a filling factor N^/Ns, where Ns = Ne. + '\- 
would describe single-hole tunneling in an almost filled Landau level. The program used 
for the calculations performed here is also capable of calculations near the filling factor u = \. 
Another interesting proof for the tunneling would be to prepare the system in a state where the 
quasihole is localized at one delta potential and to evaluate its time evolution. An oscillation 
of the quasihole between the localized states should be observable. 
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6 Inhomogeneous systems: Constrictions in the FQH regime 
6.1 Modeling the constriction 

After having seen in section |3 that the quasiholes indeed show tunnehng between two delta 
potentials, the aim is now to investigate their behavior if they are not bound to point-like 
potentials but are allowed to move in a system with a constriction. We want to infer their 
transport properties in a fractional quantum Hall system with a constriction directly from 
their motion in the time domain. 

The system to be considered therefore not only subjects the electrons to a homogeneous 
magnetic field, but also to some external potential V{x,y). This potential serves to model a 
constriction for the electrons. Its shape is that of a wall that crosses the system parallel to 
the y-axis with a notch of controllable width and depth inside this wall which forms a sort of 
a passage through it. 

Two different kinds of potentials are used here. The first one is a Gaussian shaped barrier, 
the second type is a "delta" barrier in k-space. The first type is just a potential in x- 
space. Since the calculations are to be performed in the periodic basis derived in section 
I3.1..S( the potential as well has to maintain the discrete translational symmetry in order to be 
treatable in this basis. Put differently, V{x, y) has to commute with the magnetic translations 
found in equations © and (jlOf) . Since V{x,y) only contains the operators x and y, every 
periodic potential with a period of one unit cell (in x- and in y-direction) commutes with these 
translations. The part of the potential creating the "wall" parallel to the y-axis therefore has 
to be made periodic in the x-direction. This is done by a superposition of Gaussian profiles, 
each of which is located in one unit cell. 



centered around xq in the x-direction in every unit cell. The dimensionless parameters s and 
w control the height and the width of the potential respectively. The single-particle matrix 
elements for the basis functions (|2()j) in the lowest Landau-level (n = 0) can be calculated 
straight forwardly and are found in appendix El Equ. ((^. 

These matrix elements are, as already anticipated from the translational invariance in y- 
direction of Vj^jaiiix), diagonal in k. Therefore in the many particle system, the potential only 
causes a mixing of states in subspaces with constant total y-momentum. Further on, these 
elements are real, which is an advantage for the numerical diagonalization. 

To "cut" a notch into this wall, a second potential Vnotch is superposed. In x-direction it 
inherits the shape of V^aiii in y-direction it is similarly Gaussian shaped. It is defined as 




(65) 



Its height is given in the overall energy unit 



defined by the coulomb interaction and it is 




(66) 



The matrix elements can again be computed exactly and one ends up with Equ. ()94() . appendix 
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The shape of a typical potential can be seen in Fig. 1241 The colored lines are equipotential 
lines, the color code represents the height, x- and y-axes are the coordinates in the unit cell. 
One technical point is the adjustment of the depth of the notch in the barrier: For a width 
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Figure 24: Gaussian shaped potential V{x,y) = Vwaa{x) + Vnotchix,y) of a constriction ac- 
cording to Eq. 1)65(1 and ()66() . The contour depicts equipotential lines. 

larger that approximately 0.3^2 the Gaussians of neighboring cells begin to overlap. By this 
the effective peak height at yo is higher than defined by Snotch- This overlap has to be taken 
into account when setting the parameters (in our program this is done automatically) to keep 
the potential positive definite. 



6.2 Corrections from the next Landau level 

As described in section l4.1.41 the many body problem will be treated in a finite basis containing 
only states in the lowest Landau level which is equivalent to diagonalizing the projected 
Hamiltonian PlllHPlll- While the crucial properties such as the gapped ground state, 
the incompressibility, the correlations and quasihole excitations were seen to be reproducible 
in this truncated basis, a problem occurs when evaluating expectation values related to the 
kinetic momentum operators. This will be described in the following. 



6.2.1 Vanishing kinetic momenta 

In a strong magnetic field an electron is expected to respond to an applied in-plane electric 
field by a drift perpendicular to E and B. This E-cross-B-drift velocity Vd^y is caused in our 
case by an external potential V{x) = enuVo(2;//o), which is measured in the overall energy 
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unit enu = ^ defined by the Coulomb interaction. So we obtain 



1 dV 



VdA 



B dx 
dVQ 

elQUJcfn 3(x//o) 
e2 dVn 



h d{x/lo) ' 

This drift velocity thus is important if the external potential varies on our energy scale over 
length scales comparable to the magnetic length. It will be seen, that for creating effective 
tunneling barriers, the potential will be of the order of the Coulomb interaction and it varies 
from zero to its peak value within a few Iq. So the drift velocity is not negligible. 

On the other hand, calculating the expectation value of the kinetic momenta for a state 
within the lowest Landau level returns zero. To obtain non-vanishing contributions, a correc- 
tion from higher Landau-levels is necessary, since these operators are identically zero within 
the lowest Landau level. 

Vanishing kinetic momenta Hx and Uy in the lowest Landau level formally means that 
PhLL^mPLLL = 0. Equivalently, this can be expressed as vanishing matrix elements (0, A;|nj |0, j) 
where i = x,y. 

Taking the wavefunction from (|19|) and using that Tlx and Ily commute with tx{Xp) and 
making use of Equ. to express the kinetic momenta as Hx = \J^^^i{a) — a) and liy = 

^^Y^{a) + a), the statement above is appreciated as follows 

P 

{^Aprb^x/y\^,j)prb = (0, /clp^fen^/j, V exp(-i/ca) ( (27rj + /?) --^ + /cLi) |0, 0) (67) 

keZ 

P 

= (0,Z|p,b Vexp(-iA;a) t4(27rj + /3)-^ + /cLi) n^/y|0,0) 

fcGZ 



oc 



/2 

(0, l\prb V exp{-ika) tx ((27rj + (3)-^ + kLi) (at ± a)|0, 0) 

(0, Z|prb|l, j)prb = 0. 



This means that the quantum mechanical current operator only yields non-vanishing contri- 
butions if the wavefunction has contributions in higher Landau levels. This will be analyzed 
more precisely in the following sections. 

6.2.2 Perturbative contribution of the first Landau level 

To solve the contradiction between the vanishing kinetic momenta and the expected E x B 
drift, the amount of admixture of the second Landau level due to the electric field caused by 
the potential is calculated. 

This is a single particle effect and hence it is sufficient to consider the single particle 
Hamiltonian which can be separated in the kinetic part like in Q - here called Hq - plus a 
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perturbing external potential V, which will be accounted for by perturbation theory. Let V{x) 
be a potential only depending on the x-coordinate that additionally possesses a periodicity 
with respect to the unit cell's size. Thus the wavenumber k remains a good quantum number, 
meaning that there is no mixing of eigenfunctions with different k caused by V{x). The 
eigenfunctions of the unperturbed system Hq are known from ()19() . The first order correction 
in perturbation theory now reveals 



^-^ {71, k\V{x)\m,k) 



1 

flLOr 



(-(0,A:|y(x)|l,fc) |0,A:) + (l,fc|y(x)|0,A:) 



(68) 



In the second step all higher Landau- levels except for the first two were neglected. The 
appearing matrix element can be calculated using the basis in direct space from H2U() . Due 
to the diagonality in k all wavefunctions only have strong contributions in the vicinity of 
—Xk- So a Taylor expansion of V{x) around —X^ makes sense. The linear term will cause a 
transition between neighboring Landau levels. This follows from the orthonormality relation 
J^^dx exp{—x'^)Hn{x)Hm{x) = (5n,m2"n!-v/vr for Hermitian polynomials. Thus, all higher 
terms of order n produce contributions only between Landau level m and m ziz n and can 
therefore be skipped without increasing the error of this approximation. The remaining 
matrix element can be calculated analytically, where the Taylor expansion of V{x) is used. 
This leads to the matrix element 



{0,k\Vix)\l,k) 
Using that, the correction of |0,A;) becomes 

' ^ V2 dx/k 



k dV 
y/2 dx 



{-\Q,k) + \l,k)). 



(69) 



(70) 



The formula (|7U() indicates that the correction of the wavef unction however is small if V does 
not vary much in units of huJc over distances of order of the magnetic length Iq. To lowest 
order, the wavefunction remains normalized when adding the correction term. So it is possible 
to calculate the expectation value of the current in y-direction in this new state as 



((0,A:| + 




0,k) + 



{0,k\Uy\l,k) 



(71) 



In the first step the vanishing matrix elements of Hy in the same Landau-level was used as 
well as the hermiticity of this operator. The second step involves the explicit calculation of 
(0, A;|nj^|l, k) = in direct space. Thus we gain the result for the drift current which was 

classically anticipated from the E x B drift. This correction, as stated above, is important if 
V varies on the energy scale enu over a distance of order Iq. 
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6.2.3 Hellman-Feynman Theorem 

Another way to appreciate the correction from the last section is via the Hellman-Feynman 
theorem. This theorem is applicable to parameter dependent eigenvectors of a Hermitian op- 
erator depending on the same parameter. The parameter in question here is the phase factor 
(5 from the periodic boundary conditions in Equ. p4|) . which is useful when calculating the 
kinetic momentum. The parameter appeared originally in the eigenvalue equations for the 
periodic boundary condition (Equ. (1^1 . but however it is possible to apply a unitary transfor- 
mation like in section 113.1.51 Equ. (|21() . which causes (5 to reside in the Hamiltonian's kinetic 
y-momentum, as shown in Equ. (|23j) . By this transformation, the transformed eigenvectors 
obtain a fixed phase factor /?' = for the periodic boundary condition in y-direction. Due 
to its unitarity, the transformation does not affect any observable, especially all eigenvectors 
are transformed by the same transformation U{y) to become eigenvectors of the transformed 
Hamiltonian iif^ = U{y)HW (y). After transforming, the Hellman-Feynman theorem (for 
example from [201) applicable to this problem. 

Now assume we have the same system as in the previous subsection — namely a single- 
particle problem with a potential V(x) — and have calculated the eigenvectors of the Hamilto- 
nian Hp which are denoted by \k, f3). Normalized eigenvectors assumed, we have -^{k, (3\l, (i) = 

0, from which follows = — (fc, j3\ ^^^^ . In our case we yield for ff^ 

{k,m'pV,(i) = Ek,f,6k,i (72) 

dEk,p ^ , „, d\l,l3) M^^'Pu m 

-Qp-^k,i = {k,(3\—^{Ek,fi-Ei^p) + {k,l3\-^\l,/3). 

Where 5k^i denotes the Kronecker symbol. Keeping in mind the /3-dependence of H^ ^, the 

derivative of H'^ can be expressed as = ^^H^^^. Thus for k = I we arrive at a connection 
between the energy dispersion and the current as 

-l-{Uy), = dpE,,p. (73) 

What we actually calculate when restricting the basis to the lowest Landau-level are not 
the eigenvalues of H'^ but those of PlllH'^Plll where Plll is the projection to the lowest 
Landau-level (see section l4.1.4|) . Further on, because V{x) only depends on x, this potential 
does not mix any of the eigenstates in the lowest Landau-level. That's why the basis states 
are already eigenstates of the system, only their degeneracy is lifted and replaced by the 
following dispersion 

Ek,0 = ^ + {kf,\V{x)\kp) (74) 



ildv 

Lo dx 



-X,. 



k.0 



where in the second step the localization of the basis states around —X^^p as seen from Equ. 
()2U() was used together with the assumption of only small variations of the potential V{x) on 
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the length of Iq. In the last line the definition of ^fc,/? from the same equation was taken into 
account. Merging Equ. ((75|) into (ffl]) . we gain the same result as obtained in the previous 
section by means of perturbation theory. If, on the other hand, this perturbation is not taken 
into account, Hy will be zero (according to section [6.2.1(1 which will lead to a contradiction 
between Equ. 1731 and ITU In the case of the current operator, the contributions from the next 
Landau level are necessary and sufficient to arrive at a consistent picture for the projected 
system. 

6.2.4 Corrected kinetic momenta 

In the previous sections it was shown that consistency is achieved only if the drift currents 
due to admixture of the next higher Landau level is taken into account. The reason for this 
simply lies in the nature of the perturbation expansion. Restricting the wavefunction to the 
lowest Landau level means that it is the solution of the perturbation to zeroth order \0,k). 
On the other hand, the eigenenergies we calculate are eigenvalues of a Hamiltonian containing 
the external potential V{x) and especially the diagonal terms (0, fe|y|0, k), which are already 
corrections of first order in perturbation theory for the energy. Thus it is not surprising that 
inconsistency appears if we evaluate an operator {Ilx,y) in a state that was obtained by taking 
the limit hcoc ^ enu. The correct procedure is to first take into account mixing of higher 
Landau-levels by perturbation theory, calculate the expectation value and take the limit at 
the end. The reason for this is, that the operator itself can have matrix elements whose 
contribution between adjacent Landau levels diverge as we take the limit. 

In the single-particle system it is easy to correct for this lack, as shown in section l5.2.2l Since 
the external potential is the only perturbation here, instead of working with the eigenvectors 
in first order perturbation, we can directly put this correction into the single-particle operators 
of which we wish to evaluate the expectation values, namely the current operators. Reviewing 
Equ. ((7T|) . we obtain the same result for the expectation value if we correct the operator IIj^ 
with the following additional term 

|0, A:)(0, (75) 

6.2.5 Corrected current density 

The results above were only shown for the y-component of the momentum. But clearly the 
corrections cannot depend on the gauge, thus they must taken into account also for the x- 
direction. What we are interested in here is not the expectation value of the momentum but 
rather the quantum mechanical current density, which is (as a single-particle operator) 

J(ro) = ^ (G5{f - fo) + 6if - ro)G) , (76) 

with G being the velocity operator. As seen from the sections before due to the projection 
obviously G ^ ^11 as expected from the correspondence principle is not valid. Also the 
calculated correction for IIj^ in Equ. ((75|) cannot be used generally, firstly because it depends 
on the gauge and therefore cannot be applied to 11^. and secondly it was assumed that the 
potential only depends on x. A more rigorous approach for deriving the velocity operator 
just using the features of the projections of x and y onto the lowest Landau level (as used by 



n„ = n„ 



E 



1 dv 

ujf dx 
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Shankar in [^) is found in appendix IB. 21 This approach yields a different kind of current 
density operator which doesn't have any contributions from the kinetic momenta but is similar 
to the correction term calculated by perturbation theory 



Due to its agreement with the corrections calculated before and its consistency with Ehren- 
fest's and Hellmann-Feynman's theorem, this current density is believed to be the correct 
one. Another feature can be directly deduced from its formula (|77|) . The current density is 
perpendicular to the potential's gradient. By this, the electrons flow along the equipotential 
lines as expected in high magnetic fields. This will be seen in Fig. 1251 

6.3 Quasihole excitations in inhomogeneous Systems 

After having confirmed the stability of the quasihole states in section 14.2.41 for Coulomb 
interaction and in section 14.2.31 for hard core interaction and having seen the additional 
feature of being a ground state of the homogeneous system in the latter case, we now want 
to inject a quasihole into an inhomogeneous system and calculate its time evolution. First 
we are considering a weak potential that actually is no real constriction for the electrons but 
which for sure doesn't destroy the correlated state. 

6.3.1 Weak potentials 

For a system of 5 electrons and 15 flux quanta with hard-core interaction, different soft-walled 
potentials are investigated. The shape of the potential is a wall with a notch inside. The 
depth of this notch is the only parameter which is varied in what follows. The height of 
the potential was chosen such that the fractional quantum Hall state of the system is not 
destroyed. This is confirmed by comparing the chemical potential for adding an electron 
(from Tbl. ^ of 2.91-^ to the peak height of the potential (see Tbl. IIJ). The system is in 
a regime, where electrons can simply overcome this barrier. Therefore, the density hardly 
responds to the barrier potential, which can be seen as a sign of incompressibility (compare 
Fig. I25() . Rather than as a tunneling barrier, the potential should in this case be regarded as 
a source of a dispersion E{k) for the different single particle states k. This in turn causes a 
dispersion for the many particle states which are serving as a basis. This dispersion should 
cause a quasihole to move, due to its charge. In Fig. [^S] this motion can also be found 
in the electronic current density calculated according to Equ. ((77|l. The current is due to 
the electrons drift along the equipotential lines. This can be confirmed by comparing the 
current density to Fig. 1261 Although the quasiholes have positive charge a similar motion is 
to be expected for them since not only the Lorentz force of the -B-field depends on the charge 
but also the force by the constriction potential: For electrons it is repulsive, for quasiholes 
attractive. 



A quasihole is generated at the initial position (0.2, 0.0) and its time evolution is calculated. 
Due to the additional quasihole there are now 16 flux quanta inside the system. Four different 




(77) 
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Figure 25: Left: Electronic density of the ground state for the inhomogeneous 5/15-system. 

The potential is defined by the values from Tbl|3] , the depth of the notch is 
0.04enu. Due to the weak potential there is hardly a response of the density 
(incompressibility). Right: Current density for the same state. The current flows 
along the equipotential lines. 
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values for the depth of the notch are chosen. The potential landscape for those values along 
with equipotential lines are plotted in Fig. I^H 



Height of barrier (^) 


Width of barrier (Li) 


Width of hole (L2) 


0.1 


0.4 


0.3 



Table 4: Common parameters used for all barrier potentials. Only the depth of the notch is 
being varied. 



Classically we can assume that the quasihole — treated as a charged particle of charge qqh = 
|e — should follow the equipotential line it was created on. A potential without any notch 
should lead to a downward drift of the quasihole. Shallow notches in the potential cause a 
slight deformation of the equipotential lines. As long as the equipotential line the quasihole is 
starting on does not come close to the saddle point of the potential, no qualitatively different 
picture from the case without any notch is expected. In the case of an equipotential line 
crossing the saddle point of the potential, the quasihole has two possibilities: It can pass 
through the barrier as well as it can pass by. A yet stronger notch leads to equipotential 
lines connecting the left and the right part of the system. A quasihole starting on such a line 
should therefore pass from the left half to the right half of the system. The calculation of the 
density's time evolution for an initial state with a quasihole in the position Xhoie = 0.2 and 
Uhoie = 0.0 has been performed to confirm or dismiss these expectations. 

The density of an initial state is found in Fig. 1271 This state is created for each of the 
potentials separately by application of the quasihole creation operator to the ground state 
of the Ne/Ng = 5/15-system with the respective barrier. The densities however all look 
very similar to Fig. 1271 which again is due to the small impact of a shallow barrier on 
the incompressible system. Like in section I4.2.4( the time evolution of the single particle 
occupation numbers is computed. As an illustration that goes beyond occupation numbers, 
for distinct points in time the electronic densities for the systems are calculated. 

At time t = figures [271 through 1311 look similar: The quasihole can be identified with the 
quantum numbers 1 ... 5 which are less occupied than the average (|) which results in a 
density similar to the one in Fig. 1271 Fig. [27| illustrates the occupation numbers in the 
case of no notch inside the barrier. The barrier thus is translationally invariant along the 
y-axis. This means that it cannot cause any mixing of states with different k, since the 
crystal momentum in y-direction is conserved. However, the barrier causes a dispersion in 
the single particle states E(k) which makes the quasihole drift downward in parallel to the 
barrier. This drift cannot be infered from Fig. 1271 but it was confirmed by calculating the 
density as a function of time. The formal reason for this drift is the phase difference the 
expansion coefficients acquire due to the dispersion. 

If we cut a shallow notch into the barrier (compare Fig. 1261 up right), mixing of single 
particle states takes place and causes the occupation numbers to fulfill a more interesting 
evolution, as found in Fig. 1281 Starting at the same initial position as before, the quasihole 
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Figure 26: Contom' plot of the external potentials. Left up: No notch; right up: Depth of 
notch 0.02enu; left down: Depth of notch 0.04enu; right down: Depth of notch 
Omenu. 
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Figure 27: Left: Density of the initial state for t = with a quasihole at (0.2,0.0) (look- 
ing similar for all potentials). Right: Single-particle occupation probability nk{t) 
depending on time for the system without a notch. 
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quantum number k 



Figure 28: Time evolution of the single particle occupation numbers nk{t) for a quasi hole 
initially at (0.2,0.0). Depth of the notch is 0.02enu. 
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seems to move along a periodic trajectory that has its closest approach to the barrier at 
approximately 250 . . . 300h/enu and reaches its farthest distance at t = and t ~ 550h/enu. 
This is in agreement with the expectation that it should follow the slightly curved, periodic 
equipotential line it started on. Fig. EH (right) is a plot of the density at t = 500h/enu, 
where the quasihole has almost reached its starting point again, due to the periodic boundary 
conditions. This explains the periodicity in the occupation numbers. 

Another interesting finding in Fig. |2H1 (left) are the points in time, where the signature of 
the quasihole becomes clear (t = 0, t ~ 550h/enu) and those where it appears to be washed 
out (t = 300h/enu). The density for the point t = 3007i/enu where the quasihole is not 
clearly visible in the occupation numbers reveals Fig. I29L left. The cause for the appearent 
disappearence is due to an accumulation of charge directly obove the quasihole. This amount 
of charge obviously compensates locally for the quasihole's charge deficit and therefore makes 
the dip in the occupation numbers vanish. For later times (t = 500/i/enu) this accumulation 
of charge seems to decay again, (see Fig. I29| right). 




Figure 29: Left: Density at t = 300 in a system with a shallow notch (depth 0.02enu). 

The quasihole is "chased" by some charge accumulation centered at ~ (0.30, 0.83) 
Right: Density at time t = 500g^ of the same system. The charge accumulation 
is decaying again, the quasihole almost rearrived at its starting point. 



The moving quasihole seems to create some excitations of the system as a side effect. The 
fact that the charge deficit of the quasihole is almost cancelled out at t = 300^/enu by the 
charge accumulation renders the speculation temptative to assume the accumulation to be a 
quasielectron. But apart from this indication there are is no direct verification. 



A notch of depth 0.04enu leads to an equipotential line that connects the starting point of 
the quasihole with the potential's saddle point (see Fig. I26|) . The interesting point is the 
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Figure 30: Left: Time evolution of the occupation numbers for the quasihole and a depth 
of notch 0.04enu. Right: Density at time t = 380^^ for the same system: The 
quasihole is spread out along the equipotential line. 
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quasihole's behavior when arriving at the saddle point, where it has the possibihty to pass by 
the barrier as well as it can cross it. In Fig. 1301 this point is reached at about t = 200h/enu. 
Until this time the quasihole clearly moves towards the saddle point. Shortly later the plot 
suggests that the hole splits or spreads out. At i ~ 380h/enu there are two areas, one to 
the left the other to the right of the barrier, in which the occupation numbers are lowered. 
Comparison with the density in Fig. IHUl confirms the assumption that the formerly spacially 
localized quasihole is now spread out along the equipotential line. The part left to the barrier is 
moving downward the right part is moving upward. Due to the periodic boundary conditions 
it is hard to draw conclusions for later times, since several images of the quasihole, that have 
passed the borders of the unit cell, are superposed. 

Obviously, this kind of potential causes a more complicated behavior. Namely the splitting 
of the quasihole at the saddle point cannot be understood classically. One might think of 
it as a superposition of two states, the one of which having a quasihole left to the barrier, 
the other one one to the right. However, this interpretation was not checked, but could be 
investigated by projecting the system's wavefunction ^{t) to the one or the other final state. 




quantum number k i 

Figure 31: Left: Time evolution of occupation numbers for the quasihole initially at (0.2,0.0) 
for a notch of depth O.OSenu. Right: Density for t = 280/i/enu for the same system. 
The quasihole around ~ (0.65,0.80) has passed the barrier. 



Finally, in the case in which the notch is very deep the behavior resembles qualitatively that 
of the shallow notch, only rotated by 90 degrees: The quasihole travels along the equipotential 
line and arrives on the right side of the (widely opened) barrier. At t = 280h/enu in Fig. |^ 
the quasihole can be identified as the dip in the occupation numbers 9 ... 11 . Comparison 
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with the density in Fig. |^ shows the hole which for the most part is on the right of the 
barrier. Here again, to the left and to the right of the quasihole the density shows two areas 
of charge accumulation, similar to case of the shallow notch mentioned earlier. 

To end with, the following conclusions can be drawn from this section. Shallow potentials 
can be regarded as sources of dispersion, but do not affect the density of the incompressible 
state to a big extent. This dispersion leads to an E' x S-drift of charged quasiholes. Inserting 
a quasihole on an equipotential line that doesn't come close to a saddle point of the potential, 
the motion of the quasihole coincides with the classically expected E x i?-drift. 

What exactly happens at the crossing point of equipotential lines is unclear, but a superpo- 
sition of two states both carrying a quasihole on different sides of the barrier seems plausible. 
As seen from the system of two delta potentials, tunneling of quasiholes is possible. The 
situation found here is not much different, since classically the quasihole must strictly stay 
on the equipotential line it started on. If however equipotential lines left and right of the 
barrier come close to each other near a saddle point of the potential, communication between 
quasihole states on the left and those on the right takes place and the quasihole can partly 
pass over to the other side. Following this process, the part of the quasihole resting on the 
left side travels on in its original direction, while the part on the right moves in the opposite 
direction along the equipotential line. The quasihole being spread out along the equipotential 
line (Fig. I3UI density) resembles the behavior in the system of two delta potentials, where 
there was a high occupation probability for the quasihole between the two peaks. Further 
investigations should focus on the state obtained after the quasihole crossed the saddle point 
of the potential. This could be done by projecting it to an assumed final state. 

In these inhomogeneous systems the quasiholes suffer from dispersion effects due to their 
spacial extent, which causes a deformation of the quasihole. Nevertheless, in systems where 
the quasihole does not come close to any saddle point of the potential, it did not decay 
during the time to cross the unit cell. Larger systems with shallow potentials should lead 
to less distortion, because the extent of the quasihole will be realtively smaller. A relatively 
narrower spacial extent also means that the curvature of the dispersion E{k) has a lower effect 
on the quasihole. This is because the E x B-drih of a single-electron state is proportional 
to H (section IfiT^l Eq. EH) and in a bigger system the extent of a quasihole in fc-space is 
relatively narrower as well. Thus all single-particle states "carrying" the quasihole move with 
approximately the same drift velocity and hence the deformation of the quasihole is low. 

It would also be interesting to have a closer look at the described effect of charge density 
accumulation along the trace the quasihole moves. The question here is to check if the maxima 
in the electronic density can be identified with quasielectrons or if they are charge density 
waves excited by the charged quasihole moving through the system. They were found to be 
created spontaneously and also decayed much faster than the quasihole (which is stable, apart 
from its deformation). A quasielectron therefore seems to be the less plausible explanation, 
since we would expect its features to be more similar to those of a quasihole. 

6.3.2 Strong potentials 

The weak potential in the previous section did not force the density of the system to zero 
inside the barrier. Therefore it is not an appropriate choice to model a constriction like a 
point contact. To create an effective barrier for the electrons a stronger potential of the order 
of their chemical potential is needed. 

A system with a barrier enforcing a strong dip in the density is investigated here. This 
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barrier can not yet be regarded as a tunneling barrier since the density below it is still too 
high. Compared to the chemical potential of 2.91^, its height is roughly |. The u = ^ 
filled system again has 5 electrons and 15 flux quanta. The electrons interact via hard-core 
interaction. An impression of the potential can be gained from its contour plot in Fig. 1321 



Height of barrier (^) 


Width of barrier (Li) 


Depth of hole (g) 


Width of notch (Lg) 


0.5 


0.2 


0.25 


0.05 



Table 5: Parameters of the potential to create a tunneling barrier. The parameters are defined 
in Eq. (l65|) and (IHHl). 
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Figure 32: Left: Plot of the potential with parameters given in Tbl. [51 Right: Resulting 
electronic density of the ground state of the 5/15-system with hard-core interaction. 

The eigenenergies of the ground states of a system with 15 and 16 flux quanta (see Tbl. 
respectively, show that in a system with a barrier it is more favorable to have more flux quanta. 
This leads to a negative energy per flux quantum (adding a flux quantum lowers the energy) . 
It is understandable, if one keeps in mind that the incompressible state of a 5/15-system needs 
all the available zeros of the wavefunction to reduce the hard-core interaction and thus does 
not have any freedom to react to the barrier. The additional flux quantum in turn allows 
to flx a zero near the barrier's peak to reduce the potential energy and still has the freedom 
to reduce the interaction to zero. In the 5/16-system the expectation value {Vhard-core)GS of 
the interaction potential decreases tremendously by a factor of approximately 9 (see Tbl. 
compared to the 5/15-system. So there is a tradeoff between reduction of potential energy 
and interaction energy. 

When looking at the excitation energy of a quasihole, these energies are always given with 
respect to the ground state energy of the 5/16-system with barrier. Another point to make 
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is that, in contrast to the homogeneous system, the energy of a quasihole state depends on 
the position where the hole is created. This wih be surveyed in the following section. 



System 


GS energy (|^) 


interaction energy (l^short-range)^^ (i^) 


5/15 H 
5/16 H 


0.803164 
0.676082 


0.179325 
0.020571 



Table 6: Ground state energies and interaction energy {ysYiQ^t-Yange^GS of systems with 5 
electrons and 15 or 16 flux quanta, respectively. The potential is defined by the 
parameters from Tbl. El By inserting an additional flux quantum, the interaction 
energy can be reduced by approximately a factor of 9. 



6.3.3 Excitation spectrum in dependence of the quasihole's position 
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Figure 33: Excitation energy E for quasiholes inserted at different positions {xhoieiUhoie) in 
the unit cell. The quasihole was created by application of the creation operator to 
the 5/15-system's ground state. The excitation energy is calculated according to 
equation (|6T|) . 



Fig. Oni shows the dependence of the excitation energy on the position the quasihole was 
created. It reveals that the energy of a quasihole inside the barrier {xhoie = ^-^^Uhoie = 0.0) 
is close to the groundstate energy of the system. Holes created further away have higher 
energies. This implies, that the ground state of the 5/16-system should be a state, in which 
the quasihole is somehow distributed inside the barrier. The variation of the energy with 
the y-coordinate of the quasihole is quite low. The reason is the relatively narrow notch of 
0.05-L2 — 0.5/o- The quasihole, having a characteristic size of cannot resolve this short 
ranged potential variation. The structure in the density varying on the scale of Iq in Fig. 
corroborates this. 
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The density in Fig. l32l shows that although the barrier does not yet force the density to zero, 
it has a strong impact on the areas besides the peak of the potentiaL Here, the electrons tend 
to reside in two stripes parallel to the barrier with a highest probability close to the notch. 
The local filling factor is about 20 to 100 per cent higher than |. Whether the correlations 
— which are crucial for the fractional quantum Hall effect — still exist in the inhomogeneous 
system, can be answered by looking at the interaction energy of the hard-core interaction for 
this state. According to Tbl. IHl the interaction energy of 0.179enu takes round about 22 per 
cent of the ground state energy. The comparison to the system with 16 flux quanta shows, 
that here the interaction is only 0.021enu, which takes up only 3 per cent of the total energy. 
This again fortifies the assumption, that the system with one fiux quantum in addition has 
more freedom to react on the barrier than the system at u = ^ filling. This agrees with the 
comparison of the densities in Fig. |^ (left) of a system with 16 fiux quanta to the one in 
Fig. 122 in the case of 15 fiux quanta: The reaction on the barrier potential is much more 
concentrated near the actual potential in the first case, while in the latter case the potential 
has a more long-range impact on the density. In addition, the density below the barrier takes 
lower absolute values in the case of an additional fiux quantum. In conclusion it can be stated 
that an additional fiux quantum leads to a relaxation of the competition between interaction 
and potential energy. 

6.3.4 Quasihole excitation by a delta potential 

Applying the method to create a quasihole by means of a delta potential introduced into 
the inhomogeneous system has the same effect as in the homogeneous case: One zero of the 
wavefunction is pinned at the delta's position. Similar to the case of homogeneous systems 
with Coulomb interaction, here the quasihole state created by this method results in a lower 
energy. As an example, a hole was created at {x^oie = ^-O^Uhoie = 0.5). Its excitation energy 
above the ground state is 0.101099enti. In contrast the energy obtained for the operator- 
generated hole as shown in Fig. 1331 yellow plot at x = 0, is about O.lSenu above ground state. 
Again we find the delta potential method resulting in a lower excitation energy compared to 
the operator generated state. The reason for this is the quasihole creation operator that 
was deduced from the trial-wavefunctions for the homogeneous system. Applying it to the 
ground state of an inhomogeneous system results in a more or less good approximation for the 
quasihole state. In contrast, diagonalizing with an additional fiux quantum and the constraint 
of a fixed zero in the wavefunction allows the system to benefit from its degees of freedom 
optimally. Fig. 1341 (right) shows the density for the state obtained in this way. In the vicinity 
of the barrier it almost looks identical to the state without a quasihole (left). One difference 
is the absolute value of the density below the barrier: The density rose in the case of the 
injected quasihole. This is in agreement with the explanation, that to minimize the potential 
energy, quasiholes are pinned below the barrier. Fixing one away from the potential must 
increase the density below the barrier. 

6.4 Overcoming the incompressibility at the constriction 

The incompressible nature of the ground state in the fractional quantum Hall regime has 
an impact on the system's reaction to an introduced barrier potential. By comparing the 
density profile of ground states of systems with a barrier, one of which having a filling factor 
of u = h the other owning an additional fiux quantum, we can see a different behavior. It 
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Figure 34: Left: Density of the ground state of a 5/16-system with the potential defined by 
Tbl. [SI Right: a quasihole excitation is created by a delta potential at (0.0,0.5) in 
this system. The Excitation energy is ii^ — Eqs = 0.101099enu. The profile of the 
densities near the barrier are the same, only the absolute depth in the right plot 
is less, because now one quasihole is "missing" . 

can be traced back to the incompressibility at | filling where the chemical potential of a 
homogeneous system has a step. It is the reason for a finite amount of energy needed to 
create a pair of quasiparticles (quasihole and -electron) (see section r4.1.5() . Introducing an 
external potential into this system, these particles are created and allow the system to adapt 
(if the potential energy is higher than the creation energy of a quasiparticle pair). But, due 
to the finite energy to pay, the system is rigid. Moving off the fraction of ^ by increasing the 
number of flux quanta, the excessive quasiholes are available in the system "for free". They 
make the system become more adaptive to the external potential. The densities for such two 
systems were already investigated in section EHl Figure 021 depicts the case where the filling is 
|, whereas the left plot in figure 15^ shows the same system with an additional flux quantum. 
In the latter system, the area in which the density is lowered by the barrier is much more 
localized around the barrier potential. The system at ^ in contrast shows a more long range 
impact of the barrier. 

Another way to explain these different reactions of the systems can be attributed to the 
concurrence between the interaction energy and potential energy of the electrons. Thinking 
in terms of the wavefunction, we can quantify the interaction in case of a hard-core pseudopo- 
tential by the positions of the wavefunction's zeros in the relative electron coordinates. We 
have a vanishing hard-core interaction energy whenever it is possible to have a threefold zero 
in the wavefunction upon approach of two electrons. Clearly, an additional potential tends 
to attract zeros of the wavefunction in order to minimize the potential energy. In the case of 
a delta-peak potential it was shown that a single zero is pinned at the position of the delta 
peak. In turn, there is one zero missing to diminish the interaction to zero. But, on the other 
hand, if there is one flux quantum above | inside the system, this will be pinned by the delta 
inhomogeneity, leaving 3 zeros per relative coordinate to minimize the interaction to zero. 
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Figure 35: Dependence of ground state energy, interaction energy {Vhard-core)GS gap 
on the number of additional flux quanta (Ng — 15) for a system of 5 electrons 
interacting via hard-core interaction. 



In the following three subsections different approaches are tested for their applicability 
to create an effective constriction while maintaining the correlations of the incompressible 
fractional Quantum Hall state. 

6.4.1 Pinning flux quanta by delta potentials 

The idea arising from the previous observations is to build up a constriction out of several 
Dirac delta functions, each of which pinning one additional flux quantum. The density in- 
evitably goes to zero at the positions of the delta potentials. By placing the delta potential 
adequately, it is possible to construct a wall-shaped barrier with a notch in it. Since there 
are 3 unbound flux quanta per electron remaining, the system can still exhibit Laughlin-like 
correlations with a vanishing interaction. 

This was done for a system of 5 electrons interacting via the hard-core pseudopotential. 
Three delta potentials were located at (0.5, 0.0), (0.5, 0.25) and (0.5, 0.75) (in units of cell size). 
The number of flux quanta is varied from 15 corresponding to a filling of u = ^ in case of the 
homogeneous system till 18. If the idea is correct, in the case of 18 flux quanta the system 
will show properties of a fractional ^-state. To verify this, the eigenenergy of the groundstate, 
the interaction energy Eint = {Vhard-core)GS and the gap is plotted against the number of 
additional flux quanta in Fig. As expected, the ground state energy is reduced with 

increasing number of flux quanta. The change in energy is mostly carried by the interaction 
part. So the delta potentials enforce zeros in the wavefunction at their respective position 
on costs of the interaction. In the case of 18 flux quanta the eigenenergy vanishes again, 
indicating Laughlin-like correlations that minimize the interaction. Also the degeneracy of 
the ground state becomes threefold, as in the case of a filling of i. This was already expected 
from section 14.2.21 for the case of one delta potential and one additional flux quantum and 
also holds in this case. The gap between the ground state and the first excited state reaches 
its maximum in the case of 18 flux quanta, although this increase is not monotonically. This 
may be due to an odd/even symmetry in the number of flux quanta. The absolute value of 
the gap is still only approximately 25 per cent of that in a homogeneous system, which is 
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height of wall (enu) 


width of wall (Li) 


depth of notch (enu) 


width of notch (L2) 


1.0 


0.15 


-1.0 


0.3 



Table 7: Parameters of the Gaussian barrier used in the 5/18-system. 



0.800236enu(A^s = 1) (see Tbl. 0. (enu(iVs = 1) denoting the energy unit in the case of 1 
flux quantum). In figure OHl the densities of the ground states for the four different systems 
are shown. The top left plot depicts the case of 15 flux quanta, top right contains 16, bottom 
left 17 and finally bottom right has 18 flux quanta. A qualitative difference seems to exist 
between the systems with an even number of flux quanta (right) and an odd number (left). 
In case of odd numbers the density inside the notch is a bit lower. The density of the system 
with three additional flux quanta (bottom right) shows the most homogeneous distribution in 
the region away from the delta peaks, while the cases with a lower amount of flux quanta show 
some stripes parallel to the "barrier". Also, in the case of three additional flux quanta, the 
density reaches the value of approximately 6 in the area far away from the potentials, which 
equals a local filling factor of | again. Finally the two-electron correlation functions defined 
by equation H78|l are calculated for one electron being positiond sitting ati?i = (0.0,0.5), 

9{Ri,R2) = - R,)6{r, - i?^)!^). (78) 

Figure |27| shows the results. As expected from the densities, again the system with three 
additional fiux quanta shows the most liquid-like correlations. This correlation function is 
essentially structureless in the areas away from the barrier, whereas the systems with a fewer 
amount of flux quanta still show some stripe-like structures well off the barrier. 

6.4.2 Gaussian constrictions and additional flux quanta 

The approach in the last section showed, that a number of additional fiux quanta enables 
the system to react to the potential and to keep the Laughlin-like correlations. A drawback 
of this method was that the density does not vanish completely below the barrier, since the 
zeros are fixed at certain positions. Another point is, that an additional quasihole in such 
a system will again lie completely inside the ground state sector and thus does not show 
any non-trivial time evolution. The delta potentials that create the barrier obviously don't 
produce any dispersion. 

Here we will combine some additional number of fiux quanta with a potential to create the 
barrier. As seen from the previous section, three additional flux quanta seem to be enough to 
enable the system to exhibit a density profile that refiects the barriers potential. Therefore 
we will stick here to fillings of 5 electrons in 18 fiux quanta. 

The barrier used here has the parameters from Tbl. [Tj The equipotential lines of this 
potential can be found in Fig. IHHl (left plot) along with the electronic density of the system 
(right plot). The density is sufficiently low below the barrier and also quite homogeneous 
away from the barrier. The purple colored isodensity lines {p = 6) indicate a local filling of |. 
Although the notch inside the barrier is quite large, the electrons hardly populate this area, 
instead they tend to localize more strongly left and right of the notch. Similar effects were 
found in the work of Krause-Kyora [3^1 for large notches. The ground state of this system 
is non-degenerate, but looking at the spectrum in Fig. 1221 shows three states that are close 
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Figure 36: Density for 5 electrons interacting via hard-core interaction in a system with three 
delta potentials at (0.5,0.0), (0.5,0.25) and (0.5,0.75). The number of flux quanta 
is increased: Top left: Ng = 15; top right: Ng = 16; bottom left: Ng = 17; bottom 
right: N, = 18. 



6 Inhomogeneous systems: Constrictions in the FQH regime 



73 




02 4 6 8 1 0.2 0.4 0.6 6 1 



itffc. 

Figure 37: Two-point correlation for an electron at Ri = (0.0, 0.5) in a system of 5 elec- 
trons interacting via hard-core interaction and three delta potentials at (0.5,0.0), 
(0.5,0.25) and (0.5,0.75). The number of flux quanta is increased: Top left: 
Ns = 15; top right: Ng = 16; bottom left: Ns = 17; bottom right: Ns = 18. 
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Figure 38: Left: Potential used in the 5/18-system with the parameters from Tbl. d Right: 
Density of the 5/18-system with this barrier and hard-core interaction. A local 
filling of ^ amounts to a density of 6 in the units used here. The parts of the 
system away from the barrier exhibit this filling (see equidensity line 6). 
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Figure 39: Spectrum of the 5/18-system with a Gaussian barrier (parameters from Tbl. [TJ 
and hard-core interaction. The ground state is threefold quasi-degenerate and 
separated by a gap from the excitations. 



ground state (enu) 


gap (enu) 


interaction energy {Vhard~core)GS (enu) 


0.18297 


0.08552 


0.01177 (6.4 per cent of GS) 



Table 8: Ground state energy, gap (between the highest of the three quasi-deg. GS and 
the first excited state) and interaction energy for the 5/18-system with hard-core 
interaction and a potential defined by the parameters from Tbl. 



in energy and which are separated by a gap from the excitations. It resembles very much the 
spectrum of a homogeneous = | system. The interaction energy from Tbl. [H] confirms that 
the short-range interaction is effectively minimized such that it only takes up 6.4 per cent of 
the ground state energy. 

These observations corroborate that the system is in a fractional quantum Hall state. Now, 
a quasihole can be inserted. As a starting position we chose (0.32, 0.0) which lies near the 
edge of local filling of |. A starting point of (0.2, 0.0) resulted in similar behavior as described 
below. If we assumed that the quasihole should just follow the equipotential lines, we would 
expect it to appear on the other side of the barrier after some time. It was found to be true 
for shallow potentials in section 16.3.11 After creating the initial state by diagonalizing the 
5/19-system with a delta potential at (0.32, 0.0), the calculation of the time evolution for this 
initial state is performed and the density of the system is evaluated for every time step. In 
figure lini through 1411 the difference of the density of the system containing a quasihole and the 
density of the inhomogeneous 5-18-system's ground state is shown for ascending time steps. 
The first of those plots shows the quasihole at its initial position (0.32, 0.0). Along with a dip 
in the density at this position, there is a lesser occupied region right of the notch. Left to it 
there is a small accumulation of electrons. For t = 15^/enu, the quasihole moves downward 
and, as also seen from prior calculations, is smeared out a bit. The darker region on the right 
side remains stationary. Until this time both sides of the system still appear to be decoupled. 



76 



6 Inhomogeneous systems: Constrictions in the FQH regime 



At t = 30/i/enu and t = 457i/enu the dip in the density right of the barrier starts moving up. 
This can be caused by a part of the quasihole following the equipotential line to the right 
half of the system (compare the potential in Fig. I,38j) . But at the same time the quasihole's 
larger part stays on the left side and continues traveling downward. Looking at the density's 
values at this time shows that the dips on both sides became shallower. This is what we 
would expect, if there was tunneling between the two edges left and right of the barrier which 
would cause the quasihole state to become a linear combination of two quasihole states each 
of which located on one edge. This is in agreement with the results from sectional Beginning 
at t = 45^/enu and being more pronounced at t = 60/i/enu the density starts becoming more 
inhomogeneous with several hills and valleys. This effect was already seen in section 16.3.11 
for weak potentials. At this time, in the lower left corner a dip appears which looks again 
like a quasihole and is more clearly visible for later times. For t = 75?i/enu the quasihole 
on the left side reaches its initial position again. This is the time, where periodic boundary 
conditions will definitely have a big impact on the results. Simultaneously, the stripe of lower 
density which was found for earlier times right to the barrier concentrates in one point near 
the notch and reaches its maximum depth around t = 90/i/enu. One very striking discovery 
was made by comparison of the density for t = 607i/enu with that of t = 907i/enu. These 
plots are approximately the inverse of each other. As said before this may be another effect 
caused by the periodic boundary conditions. 

To end this section with, some concluding remarks on the observed effects will be given. 
Developed from the previous sections, the approach to overcome the incompressibility by 
inserting additional flux quanta into an inhomogeneous system appears to be a suitable way 
to create an effective tunneling barrier in the fractional quantum Hall regime. The ground 
state of this system shows features (gap, threefold degeneracy, low interaction energy) of the 
homogeneous = 5 system and its density exhibits a local filling near ^ in about 40 per cent 
of the system's area well off the barrier. So the insertion of a quasihole makes sense since the 
system is in a fractional quantum Hall state. 

Although the density in Fig. 051 pretends an effective separation of both parts of the system, 
the time evolution shows a strong coupling of the left and the right part of the system. The 
periodic boundary condition in x-direction of course provides a connection between left and 
right but reviewing the densities in the time evolution they do not show any dynamics around 
X = 0. This is because the quasihole is attracted by the barrier and the important effects 
take place in the vicinity of the notch. For sure at t = 30?i/enu there is a communication 
between the left and the right edge causing oscillations on the right side for later times. The 
density distribution for this time looks quite symmetric and resembles a superposition of two 
quasihole states, each having the hole on one side of the barrier. As found in the tunneling 
setup with two delta potentials, tunneling of quasiholes is possible. The symmetric state that 
forms at t = 307i/enu could be due to tunneling of the quasihole between the two edges. Also 
the structures appearing on the right half of the system at t = 90?i/enu look quasihole- like 
and make quasihole tunneling plausible but have to be treated with some care, since at that 
time the periodic boundary conditions might already have a big impact. 

An interesting feature of the quasihole is that it first smears out during its time evolution, 
but for later times it tends to reshape. Since this smearing-out is most strongly pronounced 
at the time the quasihole is near the saddle point of the potential, it could be understood by 
thinking of the quasihole to distribute between the two edges due to tunneling. This is also 
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Figure 40: Time evolution of the density of a quasihole injected at t = into the inhomoge- 
neous 5/18-system. The parameters of the potential are given by Tbl. [7| From 
top left to bottom right: t = 0, 15, 30, 45, 60, 75/i/enu. The density of the inhomo- 
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Figure 41: Density of the quasihole state for t = 90?i/enu. 

supported by comparing the densities for t = and t = 757i/enu, the quasihole on the left side 
is located at its starting position (0.32, 0.0) (or at least near to it). For t = 757i/enu the right 
side of the system exhibits a dip in the density near the notch that was absent for t = 0. This 
is another indication that ai t = 75h/enu the system is in a superposition of states in one of 
which the quasihole is right of the barrier. This it what would be expected for tunneling. 

Not only the quasihole itself undergoes a process of decay and reformation but it is accom- 
panied by the appearance of charge accumulations which were already found in the system 
with a weak potential (section lfi.H.lj) . The startling resemblance of the density for t = 60?i/enu 
with the inverted density for t = 90?i/enu suggests that the quasiholes and these charge accu- 
mulations should be treated equivalently. To check if these accumulations are quasielectrons 
of the system they could be compared against the trial wavefunction given by Laughlin 
Since these excitations appeared in systems of moving quasiholes only, the simplest case of a 
quasihole driven by a homogeneous electric field (which can be realized according to section 
I.S.2jl in a homogeneous system should be investigated. 

6.4.3 Delta constriction in k-space 

Already in the work of Krause-Kyora [SHI it was seen that Gaussian shaped barriers tend 
to cause vast areas of depletion in the vicinity of the constriction and only small areas of 
the system remain homogeneous. This observation led to the request of finding an effective 
constriction that forces the electronic density to zero while having only small impact on 
the "bulk" area of the system far away from it. In the case of Gaussian barriers we tried 
to compensate for this by additional flux quanta which are trapped by the potential and 
were shown to help maintaining the correlations. Increasing the number of flux quanta also 
increases the basis' size and thus puts a stronger limit on the number of electrons we can 
treat. So it would be desirable to have a method which is less "flux-consumptive". 

The ansatz used here is defining the constriction directly by its matrix elements. Due to our 
gauge it is easily possible to create a localized wall by this. The matrix elements of this wall are 
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such that they hft the energy of one selected single-particle state |0, jo) in the lowest Landau 
level, whereas the energies of the others remain the same. This "potential" is actually the 

2 

projector onto this one lifted state, Vs = s^-^^-IO, jo)(0, jo|. Hence it is rather a pseudopotential 
as there is no local representation in x-space since the projector is non-local: There are two 
coordinates to integrate over to calculate its matrix elements. The diagonal terms of this 
"nonlocal potential" are given by a Gaussian shaped potential centered around x ~ ~^jo 
with Xjg = + 2ttN • '^^^ intention of this procedure is to create the narrowest possible 
barrier. Due to the localization of the single particle states on a range of the cyclotron radius 
Iq, this gives a lower limit for structures' sizes the system can resolve. Thus even narrower 
barriers than this "delta" -barrier cannot be constructed. Its matrix elements are obviously 



To create a notch inside this barrier, we take the matrix elements of a potential wall parallel to 
the y-direction as in equation (|8fl|) and multiply them symmetrically with the matrix elements 
from equation ((7^ 



Its width and depth is defined by the width and (negative) strength of this wall. The resulting 
matrix elements are found in appendix^ equation (|96|) . A system with 5 electrons and 16 
flux quanta was diagonalized using a delta barrier as a constriction. Its height was chosen 
to be lenu. In figure the density of the system is depicted. The area between x = 0.4Li 
and X = O.6L1 is depleted which is a width of 21q. In the area x > O.8L1 and x < O.6L1 
which makes up 40 per cent of the whole system the density is rather homogeneous and at 
a local filling factor of |. The two point correlation function for one electron fixed at the 
point of largest distance from the barrier is shown in the right plot of this figure, too. The 
correlation hole around the electron is well established and the other electrons are most likely 
to be found in the region approximately 21q and further away from the first one. Apart 
from the dip at the barriers position this correlation function is structureless — as it has to 
be for a liquid-like state. The ground state energy is 0.00311enu of which are 0.00204enu 
interaction energy which is quite low. These findings make this type of barrier appear to 
be another adequate system for inserting quasiholes. From a point of view concerning the 
numerics, compared to the Gaussian barrier in the previous section where three additional 
flux quanta were needed, one additional flux quantum (= one additional single-electron state) 
is sufficient here to obtain similar features, which of course reduces the basis' size and makes 
larger systems possible to be calculated. 

6.5 Conclusions on the inhomogeneous system 

The inhomogeneity induced by an external potential demands a careful treatment in taking 
the limit of high magnetic field. Specifically, corrections to the current operators are crucial 
to obtain a consistent treatment of the system within the lowest Landau level (section l(i.2j) . 

The first part of section 16.31 dealt with weak inhomogeneities, where the incompressibility 
of a = I system counteracts the response of the system to the imposed potential. Such 
a potential can be regarded as a source of dispersion causing the E x B drift of quasiholes 
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Figure 42: Left: Electronic density, right: 2 point correlation function for a 5/16 system with 
a delta barrier 

injected into the system. The time evolution of a quasihole state allows for the following 
conclusions: 

• Far away from the saddle point of the potential the quasihole behaves like a classical 
particle in a strong i?-field. In the electric field generated by the potential it shows the 
E X B drift. Although it spreads out a bit while it moves along an equipotential line it 
remains identifiable as a dip in the electronic density. 

• Equipotential lines with opposite direction of propagation meet at the saddle point of 
the potential. At this point the quasihole is smeared out strongly. It can be interpreted 
as a superposition of a quasihole moving along the original equipotential line and one 
having passed to the other side of the potential's peak moving in opposite direction. 

• While the quasihole moves through the system it induces excitations that become appar- 
ent as accumulations of electrons. The amount of charge accumulated just compensates 
the charge of a quasihole. This suggests the question whether these excitations could 
be quasielectrons; however their short lifetime (compared to the quasihole) rather con- 
tradicts to this interpretation. 

On the other hand strong potentials, that affected the density of the system and created 
an effective barrier for electrons, destroyed the correlations of the homogeneous = | state 
that are essential for quasiholes. Much of the subsequent work was dedicated to finding a 
proper means of introducing a barrier into the system without destroying the correlations. 
Doing so, some results of Krause-Kyora [HHl were helpful. The position-dependent excitation 
energy of quasiholes confirmed these particles to be attracted by the barrier. Additional 
flux quanta were found to be bound inside the barrier. A comparison of properties of an 
inhomogeneous system at the filling factor i/ = | to the case of having an amount of excess 
flux quanta, revealed the latter system to be more desirable: Its ground state shows more 
signs of a fractional quantum Hall ground state than the one of the system at fractional 
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filling. In conclusion one can state that the excess flux quanta relax the competition between 
correlations (interaction energy) and potential energy. 

In section 16.4.11 the previous observations were developed further by constructing a barrier 
through pinning of excess flux quanta at several delta potentials. Although the properties 
(gap, threefold degeneracy, liquid-like correlation, Egs = 0) of this system's ground state 
reproduced the features of the homogeneous system's u = ^ state, it is not believed to be 
suitable for investigating quasihole tunneling because of the "pathologic" barrier: It does not 
cause any drift for the charged quasiholes. 

Section [6. 4. 21 uses the best realization of a tunneling barrier found so far: A combination of 
a Gaussian barrier in a system with a relaxed incompressibility due to additional flux quanta 
leads to a ground state of the system that resembles very much a fractional quantum Hall 
state. The time evolution of a quasihole injected into this system could be studied. Although 
the quasihole's deformation was found to be more severe than for shallow potentials, it was 
still possible to identify it as a dip in the electronic density most of the time. Near the meeting 
point of the edges the quasihole — initially traveling downwards in the left half — favored 
a symmetric state between the left and the right edge. For later times the impact on the 
right half of the system evolved into a quasihole-like dip in the density. A tunneling process 
between the edges seems to be a plausible explanation for this behavior. To corroborate 
this, it would be possible to project the state of the time evolution onto a flnal state with a 
quasihole on the right side. 

Excitation of charge accumulations moving through the system, as observed earlier in the 
case of shallow potentials, were also found here. A striking symmetry between the quasiholes 
and these "hills" in the density was observed but cannot be understood so far. To check 
whether these excitations are quasielectrons would be possible by comparing them to the 
quasielectron wavefunction of Laughlin [21 • Charge density waves excited by the moving 
charged quasihole may be an alternative explanation. To clarify this, a homogeneous system 
with a quasihole driven by an electric in-plane field (realizable according to section would 
be a system to focus on. 

Finally in section 16.4.31 an alternative approach to realize an effective tunneling barrier in 
the fractional quantum Hall regime by a delta-barrier in fc-space is presented. It has promising 
features and tunneling of quasiparticles could be investigated in this system in a future work. 
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7 Conclusions and perspectives 

The motivation and the aim of this work was to consider single-quasiparticle tunnehng at a 
quantum point contact (QPC) in the fractional quantum Hall regime. Focussing on a finite 
electronic system the behavior of inserted quasiholes near this constriction should be studied 
by treating the many-particle Hamiltonian by exact numerical diagonalization. 

Initially preparatory work was necessary to create a suitable system in which quasihole 
tunneling would be possible to observe. The preparations targeted in two directions: Firstly, 
a proper means of describing an injected quasihole in our system had to be found and its 
validity checked. Secondly, a system with a quantum point contact had to be realized. For 
the second problem it was possible to resort to results of a preceding work by Krause-Kyora 

m- 

Finally the combination of these efforts allowed the main question of quasihole tunneling 
through a QPC to be addressed and results corroborating the tunneling process were found. 

In order to relate and to compare the numerical approach to known results for homogeneous 
systems, in section ITT] Laughlin's trial wavefunctions were introduced and their uniqueness in 
a system with a short-ranged interaction was pointed out. Handling a homogeneous system 
numerically in section 14.1.51 revealed the expected gapped ground state and its incompress- 
ibility for both kinds of interactions - Coulomb and short-ranged. Confirming Laughlin's 
correlations in section 14.1.61 for our system with short-range interaction showed the vanishing 
interaction energy ( ^short-range to be an indicator for these correlations. 

Addressing the injection of a quasihole in our framework, (section |4.2[) a quasihole creation 
operator was derived based on Laughlin's trial wavefunctions. Excitations it produces in 
homogeneous systems were verified to have the desired properties such as fractional charge, 
low excitation energy and stability for the short-ranged and for Coulomb interaction (section 
14.2.31 and l4. 2. 4() . An alternative approach to create quasiholes proved to be advantageous for 
Coulomb interaction, where it resulted in a lower excitation energy and improved stability. 
A side effect of this method was the finding of a localized ground state at a delta potential 
in a system with one excess fiux quantum, i.e. Ng = 3Ne + 1. 

Treating such a system with one excess fiux quantum and two delta potentials in sectional 
in an effective single-quasiparticle picture was qualitatively in agreement with the numerically 
obtained results. Section l5.ll analvzed the ground state of this system. It is predominantly 
a symmetric superposition of the quasihole's ground state at the one delta or the other, 
respectively, and complies with the simple picture of a single quasihole tunneling between 
bound states. Efforts were made in section fF).2\ and 15.31 to extract the tunneling strength t 
from the distance dependence of the tunnel-splitting. However, only the increase of t with 
decreasing separation between the potentials could be understood qualitatively. Quantitative 
statements can possibly be drawn from a similar analysis with different parameters (see 
conclusion on section EJ. As an alternative treatment of tunneling it would be interesting 
to prepare this system in an initial state with the quasihole at one delta and perform its 
time evolution. Localization of the quasihole at either potential alternating in time can be 
expected. 

Turning to inhomogeneous systems, the projection to the lowest Landau level, considered 

2 

as taking the limit hcoc ^ demands a careful treatment of evaluating the expectation 
value of the kinetic momentum operators. A consistent treatment is achieved in section 16.21 
making a redefinition of these operators necessary which is equivalent to taking into account 
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mixing with the next Landau level. 

Inhomogeneous systems were treated in two different limits: Weak potentials that leave the 
incompressible system's density nearly homogeneous and strong potentials (comparable with 
/^electron) create abrupt edges in the system and are more problematic [201 • Section 
16.3.11 deals with quasiholes in systems with shallow potentials which are easy to handle. 
Being the source of a dispersion for the single particle states, they cause charged particles 
to undergo an E x B drift. Due to their charge, quasiholes follow the classical motion along 
equipotential lines as long as they are far away from a saddle point of the potential. They 
remain relatively stable during time spans of interest. At the saddle point equipotential lines 
with different direction of propagation meet. The quasiholes' quantum mechanical behavior 
becomes apparent here and it is interpretable as a superposition of two states: One, carrying 
the quasihole along its original equipotential line, and a second where the quasihole moves on 
the corresponding equipotential line on the other side of the saddle point. 

In the limit of strong potentials, the competition between correlations and potential energy 
prohibits creating an effective tunneling barrier while simultaneously maintaining the corre- 
lations of a homogeneous system. This can be understood by thinking of the constriction to 
bind vortices which are in turn "missing" to establish Laughlin's correlations. Section [6.3.21 
corroborated this intuitive picture. Having more flux quanta than needed for realizing a filling 
factor of |, these "free" vortices are bound by the potential. In section [6.4.11 the necessary 
amount of excess flux quanta to create a QPC was estimated. 

Applying the idea of excess flux quanta to a system with a Gaussian barrier in section 
16.4.21 resulted in a usable setup to obtain abrupt edges while conserving the correlations: 
A threefold quasidegenerate gapped ground state was obtained similar to the homogeneous 
u = ^-system. Inserting a quasihole into the left edge of this system and performing the 
time evolution allows the following conclusions to be drawn: The quasihole suffers more 
severely from distortion compared to the case of a weak potential. Similar to the shallow 
potential, the motion follows the equipotential lines. When the quasihole reaches the point 
where both edges are close to each other a symmetric state between the edges is adopted. It 
is interpretable as a superposition of two states, each having a quasihole on one edge. For 
later times there is evidence for having a quasihole in both halves of the system. In order 
to support this suggestion, in a future work the final state obtained by performing the time 
evolution should be analyzed by projecting it to the supposed final state. These observations 
corroborate tunneling of the quasihole through the constriction and are thus in contradiction 
with recent theoretical results that found the tunneling process in the limit T — > to be 
forbidden. 

As a perspective for further work, an approach of creating the narrowest barrier possible 
in our system revealed promising results in section 15.4.31 This system would be suitable for 
quasihole injection as well. 
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8 Deutsche Zusammenfassung 

Die vorliegende Arbeit beschiiftigt sich mit dem Tunneln von Quasilochern im Fraktionalen- 
Quanten-Hall-Regime. Die Fragestellung erwuchs aus den theoretischen Beschreibungen des 
fraktionalen Quanten Hall Effekts und fiihrte kiirzlich experimentell auf iiberraschende Ergeb- 
nisse pUj. Kurz nach der Entdeckung des fraktionalen Quanten Hall Effekts (FQHE) identi- 
fizierte Laughlin [2] die elementaren niedrigenergetischen Anregungen des Systems als frak- 
tional geladene Quasiteilchen und Quasilocher und konnte damit die Inkompressibilitat des 
Zustands erklaren. Auch in darauf aufbauenden Entwicklungen jO] spielten diese Quasiteilchen 
eine zentrale Rolle und ihre Eigenschaften wurden in Experimenten jl3l 1181 1191 l2Uj und 
theoretischen Arbeiten [3 I14[ I16| I17j untersucht. Einige dieser Arbeiten basieren auf der 
Beschreibung des FQHE mit Hilfe von Randzustanden jlUl lllj die sich auf deren Erfolg in 
der Erklarung des ganzzahligen Quanten Hall Effekts |2lll] stiitzt. In diesem Zusammenhang 
werden Tunnelexperimente von Quasiteilchen (und -lochern) in Quanten Punkt Kontakten 
als probates Mittel zur Verifikation des Randzustands-Bildes angesehen. 

Ziel dieser Diplomarbeit ist es, das Tunneln von Quasilochern von einer anderen Seite 
her zu beleuchten. Dazu werden Quasilocher in einem endlichen System von Elektronen 
im fraktionalen Quanten Hall Regime beschrieben wobei die Frage des Quasiloch-Tunnelns 
insbesondere in der zeitlichen Entwicklung inhomogener Systeme untersucht wird. Hierfiir 
wird die Methode der exakten numerischen Diagonalisierung verwendet. Im^ Kapitel wird 
eine Einfiihrung in das Thema gegeben und die Fragestellung in den Kontext bestehender 
und aktueller experimenteller und unabhangiger theoretischer Ergebnisse gestellt. Die fiir 
die Beschreibung des Systems notige Basis zusammen mit den gewahlten Randbedingun- 
gen wird im IHl Kapitel abgeleitet. Anhand homogener Systeme werden im 0] Kapitel die 
Rechnungen anhand bekannter Ergebnisse (Inkompressibilitat) verifiziert und eine kurzre- 
ichweitige Elektron-Elektron Wechselwirkung eingefiihrt. Zwei Moglichkeiten zur Erzeugung 
von Quasiloch- Anregungen, basierend auf bekannten Versuchswellenfunktionen jS], werden 
abgeleitet. Die Stabilitat der Anregungen wird sowohl fiir Coulomb- als auch fiir die kurzre- 
ichweitige Wechselwirkung untersucht. Ein System mit einem lokalisierten Zustand eines 
Quasilochs bietet eine erste Moglichkeit zur Untersuchung von Quasiloch- Tunneln. Dazu be- 
fasst sich das El Kapitel mit der in diesem Rahmen einfachsten Realisierung eines Systems, 
das Quasiloch- Tunneln zwischen zwei lokalisierten Zustanden zeigt. Die Abhangigkeit des 
Tunnelns von dem Abstand der bindenden Systeme ist qualitativ nachvollziehbar. Zur Unter- 
suchung eines Quanten Punkt Kontaktes (QPC) werden im|Hl Kapitel inhomogene Systeme 
betrachtet. Beitrage des nachsten Landauniveaus werden hier wichtig um ein kosistentes 
Bild zu erhalten; sie fiihren zu Korrekturen der Stromoperatoren. Fiir schwache Potentiale 
(verglichen mit der Anregungsliicke) zeigt die Zeitentwicklung quasiklassisch erwartete Ex B- 
Drift des Quasilochs entlang der Aquipotentiallinien, doch auch eine teilweise Transmission 
des Quasilochs zwischen Aquipotentiallinien kann beobachtet werden. Die Inkompressibilitat 
des Systems steht der Realisierung eines QPCs durch ein starkeres Potential im Weg. Eine 
Moglichkeit, dennoch ein System mit einer Tunnelbarriere in einem inkompressiblen Zus- 
tand zu erzeugen wird gefunden und die Zeitentwicklung eines Quasiloch-Zustands in diesem 
System berechnet. Es findet eine Kommunikation zwischen den beiden Randern iiber die 
Barriere hinweg statt, die sich durch Tunneln des Quasilochs erklaren lasst. Das letzte Kapi- 
tel enthalt eine zuammenfassende Betrachtung und einen Ausblick auf Fragen, die in dieser 
Arbeit auftraten, aber nicht erschopfend geklart werden konnten. 



A Derivation of a quasihole creation operator 



85 



A Derivation of a quasihole creation operator 

In section [4.2.11 the trial wavefunction (Equ. ()59|)) of a quasihole excitation in a fractional 
quantum Hall state at a filling factor ^ = ^ was given along with the trial wavefunction of 
the ground state, both for homogeneous systems. Basing on these wavefunctions, an operator 
shall be derived that, applied to the ground state wavefunction, results in the quasihole 
excited state with a quasihole inserted at a position zq inside the unit cell. This operator 
will be used as a generalization to create trial wavefunctions for arbitrary systems (i.e. with 
inhomogeneities) . The quasihole is manifest in the wavefunction as a single zero with respect 
to the electrons' coordinates Zj. Since adding a zero to the wavefunction amounts to inserting 
an additional flux quantum, the final state belongs to a system with A'^s + 1 flux quanta. This 
operator will then be applied to the many-particle basis states of the system with Ng flux 
quanta. The result can be expressed in the many-particle basis with -|- 1 flux quanta. Once 
knowing this basis transformation, the operator can act on arbitrary initial states to obtain 
trial wavefunctions for quasihole excitations. 

As seen in section [4.2. II Eg u. ()59() . the relative zero between zj and zq is caused by multi- 
plication with a factor i?i ^vr — — ^Nr^) every electron j. The additional flux quantum 
needs to be accounted for by the center of mass wavefunction -F/v^, too. It must be replaced 
by a solution -F/v^+i of the single-particle Schrodinger equation given by Equ. ()44() for A^s + 1 
flux quanta. This wavefunction can be obtained from the center of mass part of the ground 
state by applying transformations that only act on the center of mass coordinate. More con- 
crete, given a center of mass wavefunction for Ng flux-quanta it is seen from (|44() that we can 
construct a valid one for A^g -|- 1 flux quanta simply by a spatial shift of Z Z + AZ and 
a momentum shift exp(iAKZ). From the two conditions in H44() connecting the wavenum- 
ber K and the zeros Zi, of the center of mass function follows, that this can be done by 
Zy ^ Zy^ -^n + and K ^ K + for odd integers n, k. This transformation can be 
expressed as a shift of the center of mass Z and an accompanying shift of the center of mass 
momentum. Additionally, A^g in the Gaussian must be replaced by A^, -|- 1. Applying this 
transformation on the whole wavefunction will not affect its relative part. Thus the ansatz 
for the operator belonging to the insertion of a quasihole at position zq looks like 

0.0^^0) = \%) (81) 

exp(- ^ )T--^-C-^)eMi^KZ)T-^iAZ)eM \ l 

L1L2 m L1L2 

^ V ' 

-/,.eK-.)^iV,+l(^+^)cxp(-^ LIX7=^) 

where AZ = —-^n — and AK = jj^n, n, k odd integers. Here we choose n = 1, A: = 1 
because we want to construct a state that does not differ much from the state we started with. 
Therefore, the smallest possible center of mass and momentum translations are chosen. As 
indicated by the curly bracket, the latter part of the operator is intended to fix the boundary 
conditions and yield the part of the trial wavefunction l|59|) that is written below this bracket. 

It is obvious that this wavefunction has got the desired zeros at zq produced by the first 
factor in 1811 In the following calculation we will apply this operator to the many-particle 
basis states of the system with A^^ flux quanta. It can be calculated separately, but it is also 
evident from this calculation, that the obtained wavefunction with a quasihole will obey the 
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correct boundary conditions, since they will be shown to be expandable in the many-body 
basis of the N^/Ng + 1-system. 



A.l Application of Ohoie to the basis states 

The many-body wavefunctions in the finite system with Ns flux quanta are represented as 
linear combinations of Slater determinants of one-particle wavefunctions from Equ. ()2U() . 
Therefore, if one wants to insert a quasihole into such a state, given its decomposition into 
these basis states, the operator Ohoie can just be used to act on these basis states and the 
result again has to be expressed in the many particle basis states of the system with A^s-|-1 flux 
quanta. If once this basis transformation is calculated, also the transformation of arbitrary 
states expanded in this basis is known. In the following this basis transformation is calculated 
explicitly. 

Since each product in the slater determinant contains single particle wavefunctions as 
factors, a Gaussian factor can be split off the many-particle basis wavefunction, like 

^'(Zl, Z^e) = /(2l, ^TVe) exp( (82) 

This is true for arbitrary wavefunctions expressed in this basis. Applying the operator Ohoie 
to this state yields 

*/foie(2l,---,^7Ve;^o) = Ohole[zQ)^{zx,...,ZNe) (83) 

oc n-^,(.^!4^K^) 

exp(KZ )f{{zi + — + — }) exp( — 

It was used, that the center of mass translation, due to being sandwiched between the two 
Gaussian factors, only affects the coordinates of / and that a translation of the center of mass 
by ^ is equivalent to a translation of each electrons' coordinate by jpj^- A factor exp{iK ^) 
was omitted because it is just a c-number which will be absorbed by the normalization. 

Every element of the many particle basis is a sum of products of iVg single-particle wave- 
functions (dOI), each of which depends on one of the coordinates {zi, zjye}- Expanding the 
Gaussian factor's exponent allows to separate a Gaussian factor proportional to exp(— x^) 
and leads to the following form for (poj 



(t>o,j{x,y) = — =-2^exp ( ^^-2 {-x + iy) -ika\ -2 1 (84) 
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An arbitrary wavefunction expanded in the many-body-basis can be written explicitly as 

'^{zi,...,ZNe) = {zi,...,ZNe\ ^ Cj^^,„jj^Jji,...,jNe) (85) 

Jl<---<ijVe 

where {zi, ZNe\ji, ■-, jNe) = det{{fj^{zi)}k,i=i..Njexp{ ) 

and Cj-^^,,,jj^^ are the coefficients in this basis. Now it is possible to apply the operator (|81|) 
to the state in This done, one arrives at 

"^hole = Oholeizo) ^ Cj^,...,jNe\h, -jNe) (86) 

jl<...<ijVe 

i{z*^-z*Q) Li , , ^0 , AZ 

""""" ~ — - ; = l,,JVe 



E C,„...,^Jet{{M^ ^^^^ o^ |^_l)exp(K4)/..(^^ + ^ + 

jl<...<jNe ^ ^ 



exp(- 



vr(iV. + l)Efii^,l 
L1L2 



This expression has to be written as a linear combination in the many-particle basis of a 
system with A'^s + 1 flux quanta. This is possible, since the wavefunction H86() satisfies the 
periodic boundary conditions. Using the periodicity of the ^-function (for example from [22]) 

by 

T?i(n + 7r|r) = -i?i(n|r) (87) 
'!?i(u + 7rr) = — exp(— ivrr) exp(— 2zm)t?i(u|t) 

where r = i-p-^ 

it is possible to show, that the wavefunction ^hoie satisfies modified periodic boundary con- 
ditions like 

27r 

^hoie{zi, ...,Zi + ...,ZArJ = eicp{ia)^hoie{zi, ...,ZNe)eyip{-—i{Ns + l)yi)(88) 

J^2 

'^hoie{zi,...,Zi + iL2,Zi+i,...,ZNj = ex.p{ip)'^ holeizi, Zn^)- 



Stated differently, equations H88() show, that the wavefunction ^hoie is an eigenfunction of 
the magnetic translation operators for a system with + 1 flux quanta, since the appearing 
exponential in the first line is compensated by the one appearing in the magnetic translation 
from Equ. ©. 

According to |^ the odd elliptic theta function of first order is defined as 

oo 

^i(n|r) = - Yl (-l)"exp(i7r(n + -)V + z(2n+ l)n). (89) 

n=—oo 

Using that, one can express the factor i?i(vr^^^j^|i-j^) fj{z) appearing in Equ. as a 
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linear combination of basis functions, namely 



X exp — vr — ; ^ 7 H h iV,a H 



5] 5] G(j-/.-n,i,zo)expH^^— a) </>jy=+i(z) 

h=0 \ne{iV.+l)Z 



:K{j,h,zo) 



where (/)^'''''^(2;) names the h-th single-particle basisvector of the system with A'^s + 1 flux 
quanta. 

Now the expression 1)86^ can be rearranged such that ^hoie is expressed as a linear combi- 
nation of many particle basis functions. To accomplish that, one element of the matrix under 
the determinant can be expressed as a product of two matrices by inserting (|^n)) into ()86() 
and applying the identity det{A ■ B) = det[A)det[B) on the determinant. 



^hoie = f^ii,...jiVedetf'{^i^(jfc,/i„,zo)</'J:+'(^z)}M=i,...,iVel(91) 

/ii<...<hiv^ ii<...<jjVe \ m=l J 

= Y Yj C'iiv.jiVedet({K(jfc,/im,zo)}m,fc=i.„ArJ X 

hl<...<hMe jl< — <jNe 



fjhole 



xdet (<A£+'(zOW=i,...,7V.) . 



Finally, the new coefficients C^°^^ in the last equation can be calculated explicitly as a 
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matrix multiplication. G{d,j, zq) is used from ()9U() . This results in 

l^hole) = C^f:..,h^^ \h,...,hNj (92) 

With ct;!.,,^^ = 



■:jNe 
jl<---<jNe 



Ky^-'[[]^^^ = detl{K{jp,hm,Zo)}p^m=l...N, 



K{j,h,zo) = ^ G{j -h- n, j, zo) exp{-i ^ a] 

ne{Ns+l)Z 



= -iQ{xo)Y,i-'^y~^^M-ina)exp{-i^{j + ^ + nNs{Ns + m x 
X exp(2^.^^^(j - N,n + {- - h)^^ + —^^ ,)) x 

Here (^(xo) is a real function of xq- It doesn't have to be computed, since the resulting vector 
will be normalized anyway. 



B Single particle matrix elements 
B.l The constriction 

All single-particle operators are evaluated by aid of second quantization, thus it is sufficient 
to calculate the single particle matrix elements to construct the many particle operator O = 
Oi,ja\aj from it. 

A straight forward calculation yields the single-particle matrix elements of the Gaussian 
constriction potential given in Equ. 1)65^ for the basis functions from Equ. ()2U() in the first 
Landau level. As seen from the potential only depending on x the matrix elements conserve 
the momentum in y-direction. 

w sr^ ( 1(3+ (3/2t: 



(0, j|VVa/z|0,/c) = s^^Y^Qy.^\--^{^—^^—^+ra-XQlL^\bj^k, (93) 



where a = ilw"^ + 



2ttNsLi/L2 ' 



The calculation of the matrix elements for the "notch" inside the constriction, the potential 
of which is given in 1)66^ . can be carried out analytically as well. A rather lengthy calculation 
yields 



2 I — 

{^,j\Vnotch\^,k) = Snatch^ w Whole ^ Yl {A{l)B{m) + A{m + l)B{l + l)) (94) 
where A(0 = exp (-i, (//2 - xo/L, - i(^ + ^)) 

andS(?n) = exp ( 27ri{k - j + Nsl)^ - ial - {-^^^ + Tr^wl^l^){k - j + Nslf 



90 



B Single particle matrix elements 



As a remark, in the case that yo = O.5L2, which is a notch in the center of the unit cell (with 
respect to the y-direction) , and a = or a = vr, the matrix elements for the notch become 
real. This is of course a benefit for numerical calculations since only a real symmetric instead 
of a complex Hermitian matrix has to be diagonalized. 

The matrix elements for the notch inside the delta-barrier are constructed by symmetric 
multiplication of those of the delta wall and those of a Gaussian wall parallel to the x-direction 
defined by Equ. H80() . Calculating them yields 

(O'il^wallJO' ^) = ^^walL^walL^Z^^^'P " ^ + ^^')^ + '"^) (9^) 

x-p(-(-HalL + ^)O--^ + ^^0^)- 
Therefore the matrix elements of the notch turn out to be 

(O'il^deltanotchlO'^) = ^ ((O.il^wallJO' ■?o)'5fcjo + (0, Jol'^vallJO' ^)'^ijo) , (96) 

where the depth of the notch is given by the negative strength of V^^jj while ^f^^ll 'width 
defines the width of the notch. 



B.2 Kinetic momenta in the lowest Landau level 

To obtain a gauge invariant expression for the kinetic momentum operator we state that 
Ehrenfest's theorem must be fulfilled 

-{Uy) = {y) = U[H,y]). (97) 
ni n 

We use this as a definition for the operator Uy or respectively Vy 

Vy = y='-[H,y]. (98) 

It is sufficient to restrict ourselves to the single particle problem given by the Hamiltonian 

H = Hkin + Vix,y) (99) 
= -L(ul + Ul)+V{x,y), 

where -fffcm is defined like in Equ. The problem shall be solved after projecting it to 

the lowest Landau level. This projection will also be done for the x and y operator following 
|30j . The "projected versions" of these operators will appear to obey a conical commutation 
relation. 

Defining the operators kx and py by 

k^ = -ihdx + eBy (100) 

Py = —ihdy 



and comparing them with Equ. @ respectively HlOj) from section llj . 1 . 21 identifies them as the 
generators of magnetic translations, x and y can now be expressed in these new operators 
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and the kinetic momenta (found in Equ. ©)• The projection to the lowest Landau level is 
performed by Plll defined by Equ. ()51() in section 11.1.41 and results in 

(101) 



X 




y 


— n^) 


XLLL 


= PlllxPlll = 


VLLL 


= PlllvPlll = 



eB 



As shown in section l().2.1l the projection of IT onto the lowest Landau level is zero. Since 
kx and Py commute with the Hamiltonian H/^-in, kx and Py don't cause mixing of different 
Landau levels and we can omit the surrounding projectors in Equ. (|lfllj) . 

Now, xlll and ulll inherit the canonical commutation relation of kx and py 

[xlll^Vlll] = ill- (102) 

The velocity operator can be calculated according to our definition from Equ. (|^5|) . This 
yields 

Vy = -^{[Hkin,yLLL] + [V{xLLL,yLLL),yLLL]) (103) 

1 

= T[^ixLLL,yLLL),yLLL]- 



The kinetic part commutes with ulll^ the potential yields the expected correction term. It 
can be calculated using [x\j^j^,yLLL] = k x^£j^j^ and a Taylor expansion of V{x,y) in x 
from which follows 



Wi-^yly] = }Zt\^x^ 



k=0 

1 d''V{x,y) 



[xlLL,yLLL] (104) 

x=0 



oo 



^-^ kl dx^ 

k=0 



x=0 



ill ^ ^LLL 



il 



dV{x,y) 



° dx ■ 

This results in the corrected velocity operator 

1 dV{x,y) 
^ mujc dx 

An analogous treatment of the velocity in x-direction therefore results in 

. _ 1 dV{x,y) 
mcoc oy 

This definition leads to a different velocity operator than the naive approach of taking the 
kinetic momenta. It not just gives us a correction but states that the kinetic momentum 11 
doesn't cause any current in the lowest Landau level at all. 



(105) 



(106) 
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C Two-particle matrix elements 

C.l General two-particle operator depending on fi — fj only 

Most of the two particle operators in this work only depend on the difference coordinate ri — r2 
between the two particles. For those cases it is sufficient to calculate the two-particle matrix 
element of the operator exp(ig(ri — f*2)). The matrix element of the Coulomb interaction, 
the short range interaction and the correlation functions can be obtained simply from their 
Fourier transforms. Since we use generalized boundary conditions we cannot just reuse the 
results of Yoshioka from |33 which are valid only for boundary conditions of Equ. (|14j) with 

a = P = 0. The wavevector q is restricted to values lying on a lattice q = (^^s, j^t^ where 
s,t £ Z since all the operators have to be periodic in the unit cell. Calculating this operator's 
matrix elements in the lowest Landau level spanned by our basis given by Equ. (|2(J|) yields 

{ji,h\exp{iq{fi-f2))\j3,j4:) = "^exp f-hlq"^ + iq^^{j2 - Ja) + ialj (107) 

^ ^jl-j4,j3-j2+Nml ^- _■ qyL2 ; 

-'1 -J^' 27r 

Z Z 

where q G 2tt — x — 
Li L2 

and 6' is the Kronecker delta modulo Ng. To evaluate a general two particle operator O with 
the Fourier transform V{q), it can be written in second quantization which then results in 



1 ' 

^=2 ^ ^ (ii,J2|exp(ig>l-f2))|j3,j4) a]^a]^a,,a,3. (108) 

C.2 Coulomb interaction operator 

The Coulomb matrix elements follow from Equ. (|l()8j) by using the Fourier transform in two 
dimensions for periodic functions in the unit cell. The calculation of its Fourier transform is 
carried out here because it was not found explicitly in the cited papers. By taking into account 
the interaction between all the images of one electron in neighboring cells, the Coulomb 
interaction is periodic in the unit cell. Alternatively we can think of the Ng electrons residing 
on a topological, flat torus and interacting with each other via all possible ways, 

^Sli ='-T. f-kLA-lL2ey- 

Its Fourier transform clearly only has components for ^-vectors satisfying q E (|jZ, |^Z). It 
can be calculated as follows 

V^Ziiq) = ^V— f^' dx dy y ,^ , ^ i ^^exp(-igr) (110) 

^""^^^^ L,L2 e Jo Jo ^ ^ \f-kLie^-lL2ey\ ' ^ > 

1 f°° f°° exp{-iqr) 



L1L2 e y_oo J-oo 1^ 



dx I dy 



27r 
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Writing the periodic Coulomb interaction in its Fourier representation then leads to the 
expression found in 

= 'UL;~^ ^ ^exp(fQr-) (111) 



C.3 Short range interaction operator 

In the case of short range interaction the Fourier transform is according to Equ. (|42j) propor- 
tional to —q^ and we just replace ^ in Equ. (ITTTI) by llq^, thus 

Vshorti^ = y^iWof exp(igf). (112) 



C.4 Pair distribution function 

To obtain the pair distribution function used in section [4.1.61 Equ. ()57() . we need the two 
particle matrix element of 5^^^(r + ri — r2) where the argument of the delta is taken modulo 
cell size, i.e. periodic as 5P^^{r) = 5{f+ kLiCx + lL2ey). It can again be calculated by 

aid of Equ. (|l()7j) . This is done by writing this delta function in Fourier space as 

S^^^{fo + r) = - — — exp(igro) exp(igr). (113) 



C.5 Two particle correlation function 

For inhomogeneous systems the two-particle correlation function in Equ. (|78() is used. Its 
two-particle matrix elements are clearly simply a product of wavefunctions from Equ. (|2()|) . 
thus 

{h.hWx - Ri)S{f2 - R2)\h,k) = ^l,,{RiWoj^{R2) ^o,j,{Ri)cl)o,j,{R2). (114) 
Its representation in second quantization is obvious then. 
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